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Abstract 

We exploit gauge/gravity duality to study 'thermal quenches' in a plasma of the 
strongly coupled Af = 2* gauge theory. Specifically, we consider the response of an ini- 
tial thermal equilibrium state of the theory under variations of the bosonic or fermionic 
mass, to leading order in m/T < 1. When the masses are made to vary in time, novel 
new counterterms must be introduced to renormalize the boundary theory. We consider 
transitions the conformal super- Yang-Mills theory to the mass deformed gauge theory 
and also the reverse transitions. By construction, these transitions are controlled by a 
characteristic time scale T and we show how the response of the system depends on the 
ratio of this time scale to the thermal time scale 1/T. The response shows interesting 
scaling behaviour both in the limit of fast quenches with TT <C 1 and slow quenches 
with TT ^> 1. In the limit that TT — > oo, we observe the expected adiabatic response. 
For fast quenches, the relaxation to the final equilibrium is controlled by the lowest 
quasinormal mode of the bulk scalar dual to the quenched operator. For slow quenches, 
the system relaxes with a (nearly) adiabatic response that is governed entirely by the 
late time profile of the mass. We describe new renormalization scheme ambiguities in 
defining gauge invariant observables for the theory with time dependant couplings. 
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1 Introduction 

Consider quantum mechanics with a Hamiltonian which depends on some external 
parameter A, 

H X = H + \5H. (1.1) 



The dynamics of the system induced by variations in A is well-understood, e.g., , see [1]. 
In particular, consider beginning with A = and preparing the system in an energy 
eigenstate of the Hamiltonian Hq. If the new coupling is turned on adiabatically, the 
system continues in an eigenstate with a time- dependent energy which simply traces 



2 



the changes in X(t). In contrast, if the coupling is abruptly turned on, e.g., A = A 0(t), 
the system would evolve forward in a complicated linear superposition of eigenstates 
of the new Hamiltonian. While the description of adiabatically evolving couplings 
is easily adapted to quantum field theory (QFT) [2], the description of the latter 
'quantum quenches' is less well understood in the context. However, it has become the 
subject of a vigorous research program motivated by the recent advances in cold atom 
experiments [3-5]. 

Gauge/gravity duality [6] provides a remarkable framework for the study of certain 
strongly coupled gauge theories. Although the most applications of this correspondence 
have been directed at analyzing the static properties of the boundary theories, there 
is no conceptual obstacle in applying this holographic framework to time dependant 
problems and in particular, to the study of quantum quenches [7]. In fact, early 
attention was given to the related question of describing 'thermalization' within this 
holographic framework [8] and motivated by connections with the strongly coupled 
quark-gluon plasma, there has been a renewed interest in this subject [9-12]. However, 
given the complexities of the bulk description of rapid changes in the boundary theory, 
numerical relativity is increasingly being applied to study these 'far from equilibrium' 
situations [13-19]. 

In this paper we begin a study of quenches in the strongly couple Af = 2* gauge 
theory [20-22] applying the techniques of numerical relativity. Recall that Af = 2* 
gauge theory is obtained as a deformation of the Af = 4 super- Yang-Mills (SYM), 
where a Af = 2 hypermultiplet acquires a mass m. For technical reasons, we will limit 
our present investigation to 'thermal quenches,' where the initial state is a thermal 
state, i.e., the Af = 2* theory is prepared in a microcanonical ensemble, and we work 
to leading order in in a high temperature expansion with m <^ T. 1 As explained 
in [25] , in such a thermal state, we can split the masses of the bosonic and the fermionic 
components of the massive hypermultiplet. Hence we investigate two separate classes 
of thermal quenches with 

Csym + A A (t) O a , (1.2) 

where we may have either the bosonic mass operator 2 or the fermionic mass operator 
O3 (with dimensions A = 2 and 3, respectively). Of course, the couplings which vary 
in time are then simply the corresponding masses, i.e., A 2 = ml and A 3 = To/. In 

1 Thermodynamics of Af = 2* plasma was discussed in [23-26]. 



3 



particular, we will consider the following specific profiles for these couplings 

Aa(*) = \ A [ 1 ± tanh(t/T) ] . (f .3) 

Note that the profile in eq. (1.3) with the + sign has Aa(£ — > — oo) = and Aa(* — > 
+oo) = A A . Hence it produces a transition from an initial thermal state in the con- 
formal SYM gauge theory to a final state of a massive gauge theory. Meanwhile the 
profile with the — sign then corresponds to the reverse transition from the massive 
theory to the conformal theory. 

We have introduced the time scale T in eq. (1.3) to control the rate of these tran- 
sitions. This will allow us to consider not only abrupt quenches (with T — > 0) but also 
transitions made over a finite time. In particular, we will also consider the adiabatic 
limit where T — > oo. We will find that the response to these different quenches will 
depend on how T compares to the thermal time scale 1/T. Hence it will be convenient 
in the following to define the dimensionless ratio a = nTi T, where T { is the temper- 
ature of the initial state. Note that while our calculations of these thermal quenches 
are perturbative in m/Tj, they are no restrictions 2 on the rate controlling ratio T{T ' . 

The holographic dual of the M = 2* gauge theory consists of Einstein gravity 
coupled to two massive scalar fields in five dimensions, as will be described in detail in 
the next section. Of course, within this framework, the initial thermal state is described 
by an asymptotically anti-de Sitter (AdS) black hole and the quenches are realized by 
varying the asymptotic of behaviour of the corresponding scalar with the given profile 
(1.3). This variation essentially excites the bulk scalar with a wave packet that falls 
down into the horizon from the asymptotic boundary. Our numerical simulations 
evaluate the evolution of these scalar excitations in detail and in particular, allow us 
to determine the response of the boundary theory, as measured by the expectation 
values of the stress energy (T^) and the operator (Oa)- Our main focus is to study 
the transitions from the massive theory to the conformal theory, however, as we will 
explain, the results for the reverse transitions readily follows from the previous case, 
in the high temperature limit. 

Before proceeding, let us summarize the results of our study: 
■ New ultraviolet divergences appear into the quantum field theory when the couplings 
are made to vary in time, as for the quenches described above. Holography provides 
a well-defined approach to regulating these new divergences and renormalizing the 



2 One caveat to this statement will discussed in section 6. 
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boundary theory, even in the presence of such time varying couplings. In particular, 
the necessary new counterterms are readily identified. 

■ The response to fast quenches with a < 1 exhibits an interesting scaling behaviour 
in a, similar to that noted previously in [9,13]. This scaling is apparent both in the 
immediate response of the observables, e.g., (0^,(t)), and in the comparison of the 
initial and final state, e.g., the overall change in the energy density. These scalings are 
particularly dramatic in the case of the dimension three operator, where the immediate 
response scales as a~ 2 ln(l/a) as a — > 0. For both operators, this scaling behaviour 
is divergent and naively suggests that an abrupt quench with a = is unphysical. 
However, our present perturbative analysis is not sufficient to reliably confirm this 
conclusion. 

■ For very slow transitions with a 1, the system essentially maintains thermody- 
namic equilibrium throughout the process and the response simply tracks the change 
in the coupling. The approach to the adiabatic limit with a — > oo also exhibits an 
interesting scaling behaviour with a. As expected the entropy production goes to zero 
in this limit confirming that the adiabatic quenches are reversible. 

■ The relaxation of the system after a thermal quench of the coupling is sensitive to 
whether the ratio of the characteristic time-scale T to the thermal scale 1/Tj is larger 
or smaller than one, i.e., whether a > 1 or a < 1. For fast quenches (with a < 1), the 
relaxation is dominated by the lowest quasinormal mode of the bulk scalar. This be- 
haviour indicates that the relaxation after fast quenches is universally controlled by the 
characteristic thermal time scale 1/Tj, as has been observed previously in [13,16,18]. 
In contrast, for slow quenches (with a>l), the system relaxes with a (nearly) adia- 
batic response and so the relaxation is governed entirely by the late time profile of the 
coupling. 

■ The holographic renormalization of the boundary theory with time-dependent cou- 
plings introduces new scheme-dependent ambiguities. This situation is rather reminis- 
cent of the ambiguities that arise in studying QFT in curved spacetime [27]. All of 
these new ambiguities only play a role while the couplings are varying and so they are 
not important after the coupling has achieved its final value. 

The remainder is organized as follows: In section 2, we review the dual gravitational 
description of strongly coupled M = 2* gauge theory, its holographic renormalization 
for arbitrary sources, and the high-temperature equilibrium thermodynamics. In sec- 
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tion 3, we derive the bulk equations of motion describing the gravitational dual of the 
mass quenches (1.3). We also compute the expectation values of the stress-energy ten- 
sor and the operator Oa throughout the quench. Our numerical procedure is outlined 
in section 4. In section 5, we provide a detailed discussion of the results of our numer- 
ical simulations. We conclude with further comments and future directions in section 
6. Further details of the analysis, as well as the description of the discretization of the 
evolution equations, are discussed in appendix A. 

2 Gravitational description of J\f = 2* gauge theory 

In this paper, we will use holography to study quenches of the masses for the M = 2* 
gauge theory at strong coupling [20-22]. However, to begin let us recall the field 
theoretic description of the mass deformation of the M = 4 SYM which yields the 
M = 2* theory. The field content of the M = 4 SYM theory includes the gauge field 
Ap, four Majorana fermions ip a and three complex scalars <f>i, where all of these fields 
are in the adjoint representation. Now the SYM theory can be deformed by adding 
two independent 'mass' terms [26] 

5C = -2 J d 4 x [m 2 b 2 + m f 3 ] (2.1) 

where 

2 = iTr(|0 1 | 2 + |0 2 | 2 -2|0 3 | 2 ) , (2.2) 

03 = -Trfifcfo- v / 2^ YM 03[0i,0 t i] + v / 2^ YM 03[0L02]+h.c.) (2.3) 

+ jjm / Tr(|0 1 | 2 + |0 2 | 2 +|0 3 | 2 ) • 

If one fixes m& = m/, the effect is to give masses to an M = 2 hypermultiplet consisting 
of Ai 5 2 and (f>i t 2, the resulting theory still preserves M = 2 supersymmetry. Note 
that beyond a fermion mass term, the dimension-3 operator contains trilinear couplings 
between the hypermultiplet scalars and 03, as well as a mass term for all three scalars. 
Note that the latter is a coupling-dependent correction induced at finite mass, i.e., 
so that rrifOs contains a contribution of order m 2 . The presence of these additional 
interactions is dictated by the supersymmetry algebra and the latter mass terms can be 
distinguished from those in 2 by their transformation properties under the 50(6) r 
symmetry of the M = 4 theory [26]. Further note that the dimension-2 operator 
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02 contains an unstable mass term for the scalar 3 , which is in the M = 2 vector 
multiplet, as well as positive masses for the two scalars in the hypermultiplet. The two 
mass contributions for 3 only cancel with rrib — mf, leaving this scalar massless in the 
supersymmetric theory. In the following study of mass quenches, we will vary m& and 
rrif independently and so we are not restricting our analysis to the supersymmetric 
theory. Our assumption will be that the structure of the two independent operators 
remains as described above. However, one should worry then that the theory may be 
unstable, 3 i.e., when m& > m/, however, we will also be working at finite temperature 
with rribj/T <C 1. In this case, 'large' positive thermal masses will be induced for all 
of the scalars, in particular 3 , which will prevent any such instabilities from arising. 

As we next describe, the holographic dual of these mass deformations is well under- 
stood [20-22]. In the dual five-dimensional supergravity, we can restrict our attention 
to a pair of scalars, which we denote a and \- This construction and its 'uplift' 
to ten dimensions was originally elucidated by Pilch and Warner [20]. In the full 
ten-dimensional type lib supergravity, the two scalars are Kaluza-Klein modes which 
deform the original AdS§ x S 5 geometry dual to the M = 4 SYM theory. According to 
the general framework of holographic RG flows [6], the asymptotic boundary behaviour 
of the supergravity scalars contains information about the couplings and expectation 
values of the dual operators in the boundary theory. Here, these are the mass param- 
eters and the two operators presented in eq. (2.1). Our nomenclature will be that the 
scalar a is dual to the 'bosonic 1 mass term 2 , in eq. (2.2), while x is dual to the 
'jermionid mass term C 3 in eq. (2.3). 

The appropriate terms in the five- dimensional supergravity action, including the 
scalars a and Xi can be written as 




de^9[\R-3(da) 2 - {d X f -V(a, X )} , (2.4) 



where the potential takes the form 



V(a, X ) = -e~ 4a - 2e 2a cosh(2 X ) + l - e 8a sinh 2 (2 X ) . (2.5) 

Implicitly here, we set have the curvature radius of the five-dimensional AdS vacuum to 
be unity, i.e., L — 1. With these conventions, the five- dimensional Newton's constant 

3 A full description of these instabilities (including their end-point) at strong coupling would require 
extending the holographic framework described below and may well require considering the full ten- 
dimensional string theory in the bulk. 
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becomes 



The resulting Einstein equations can be written as 



iV = 12cV<9„a + Ad^xduX + , (2.7) 

while the equations of motion for the scalars become 

1 9V , ldV 

Do: = -- and DX^- (2.8) 

As commented above, we are considering quenches of the mass parameters at finite 
temperature and in the regime m^j/T <C 1. This means that the dual holographic 
background will be an asymptotically AdS black hole and we will work perturbatively 
in the amplitude of the bulk scalars. In this case, the scalar sector of the effective 
action (2.4) can be expanded up to quadratic order in a and X- Furthermore, since 
the potential (2.5) does not contain a quadratic term mixing the two scalars, each of 
them can be treated independently, i.e., we study quenches in 

J ^y/=g (R + 12- i(90) 2 - \m 2 <? + 0(0 3 )) , (2.9) 



° 167rG ; 
with 

— 3 , if 6 — 2y/2\ •<=>- fermionic mass operator O s , 
m 2 = I ^ A (2.10) 

I — 4 , if = 2y/6a •<=>- bosonic mass operator 2 ■ 

The quenches which we study will be homogeneous and isotropic in the spatial 
directions on the boundary. Hence our ansatz for both the background metric and the 
bulk scalar field depend only 4 on a radial coordinate r and a time v: 

dsl = -A(v,r) dv 2 + S(v,r) 2 dx 2 + 2drdv , <f> = </>(v,r) . (2.11) 

The metric ansatz is adapted to infalling Eddington-Finkelstein (EF) coordinates, 
which has proven useful in holographic investigations of the dynamics of thermal sys- 
tems [13,14,28]. With this ansatz, we will assume that the asymptotic AdS boundary 
occurs at r — > oo with 

A(v,r)^r 2 and S^r. (2.12) 



Some limitations of this approximation are discussed in section 6. 
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Given eq. (2.11), we obtain the following equations of motion: 

2 • 39 r S - 39 r 0- m 2 

= S^S) + 2S^-2S 2 + -^-m 2 2 S 2 , (2.14) 

_L _ 



= d 2 r A- ^£<9 r £ + 4 + 0,9 r 0- ^m 2 2 , (2.15) 



= S-^ r AS + is(0) 2 , (2.16) 

= ^S + is^) 2 , (2.17) 

where we have defined for any function h(r,v), 

h = d v h+^Ad r h. (2.18) 

Note that eqs. (2.16) and (2.17) are constraint equations. In particular, the radial 
derivative of eq. (2.16) and the time derivative of eq. (2.17) are both implied by the 
first three equations. Note also that the form of the metric (2.11) is invariant under 
the residual diffeomorphism [13] 

r-»r + /(u), (2.19) 

for an arbitrary function f(v) — we use this as a check in our analysis below. 

Given eq. (2.12), it is straightforward to find the asymptotic solution to eqs. (2.13- 
2.17): 

■ When m 2 = —3, 

= Ipo + ^ (^-^oai+P^+^^-Qp'o' + ^lnr^+O^lnr), 

S = r+ l -a l + 0{r- 1 ), (2.20) 

A = r 2 + ai r + \\ - ip 2 - a[ + ^a 4 + QpoPo + t^Po - j^o) 2 ) lnr ) 
+ 0( r - 3 \nr), 

where {po,P2, ai, a 4 } are functions of the time v and we have defined h! = d v h for any 
function h{y). In addition, on the boundary, the constraint (2.16) becomes 

2 1 

-gPoa'iPo + -a'iPo a i • ( 2 - 21 ) 



Note that if this constraint is imposed on the boundary (i.e., at r — > oo) then eq. (2.16) 
will be satisfied at all values of r when the equations of motion (2.13-2.15) are satisfied. 
Finally, we observe that under the residual diffeomorphism (2.19) one has 

( Vn\ ( Po \ 

P2+Po(f 2 + fa 1 )-2fp' 
ai + 2/ 
V a 4 ) 

We emphasize that this transformation will allow us to choose a gauge where a\ — 0, 
which will do in section 3. Further, we note that the constraint (2.21) is invariant 
under the transformation (2.22). 



Po 

P2 

\ a 4 ) 



(2.22) 



When m 2 = —4, 
1 



(po+Po lnr)+C(r 3 lnr), 



E = r + -ai + C(r^ 3 ln 2 r), 



(2.23) 



A = r 2 + ai r + ^al-a[ + ^a 4 - Q^(p l ) 2 + ^p p l ^j lnr - ^(p l ) 2 In 2 
+C(r~ 3 lnr) , 

where {po,Po, Qi, ^4} are functions of v. In this case, the boundary constraint provided 
by eq. (2.16) is, 



/ 20 ... 8 . ... 10 , . 4 . 
= -2a' A - —M) + ^Po(Po) ~ ^PoPo - 9P0P0 • 



(2.24) 



Under residual diffeomorphism (2.19) one has 



/ Po^ 

Po 

dl 

V a 4 / 



/ 



\ 



Po 
ai + 2/ 

X a 4 / 



(2.25) 



Hence again, we can use this diffeomorphism to set a± — 0. Further, the constraint 
(2.24) is now trivially invariant under the transformation (2.25). 



2.1 Holographic renormalization 

The physical observables of the boundary theory are the correlation functions of the 
gauge invariant operators. In the following analysis, we will consider the simplest of 
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these, namely, the one-point functions. That is, we will examine the response of (T^) 
and ( 0&) to quenches in the mass couplings, as given by eq. (1.3). We will present our 
investigation of non-local probes of these quenches (e.g., two-point functions, Wilson 
lines and entanglement entropy) elsewhere [29]. 

The holographic computation of these one-point functions is most easily imple- 
mented using Fefferman-Graham (FG) coordinates [30], which we write here as 



dsl = % + - 2 G l3 (x k , p) dx l dx 3 , (2.26) 



and 

<P = <P(x k ,p). (2.27) 

The asymptotic AdS boundary now appears at p — > 0. As usual, to produce finite 
results in the holographic calculations, one regulates the bulk by introducing a cut-off 
surface at p = e, with some small e. However, the one-point functions extracted directly 
from the action (2.9) would still diverge as this cut-off is removed, i.e., as e — > 0, and 
thus these observables need to be holographically renormalized. The details of the 
holographic renormalization are sensitive to the conformal dimension of the particular 
boundary operator of interest, or correspondingly, to the mass m 2 of the bulk scalar <p 
in eq. (2.9). Hence we discuss the two cases for m 2 separately. 

2.1.1 Renormalization of m 2 = —3 bulk scalar 

First, we note that we are working with a flat boundary metric and thus in the FG 
metric (2.26), 

lim Gi,- = 7to . (2.28) 

Next we need to translate between the EF coordinates in eq. (2.11) and the FG coordi- 
nates in eq. (2.26). This coordinate transformation is easily constructed perturbatively 
for small p (large r) by comparing the asymptotic expansion of the bulk solutions of 
eqs. (2.13-2.15). Hence, using eq. (2.20) for m 2 = —3, we find 

v = t-p-^p 3 P0 (t) 2 + ^p 4 p (t)p' (t) + O(p 5 \np), 

r = ^(l-^P«iW+P 2 Q«'iW + ^o(t) 2 ) (2.29) 
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Note that as one approaches the boundary with p — > 0, the EF time v and the FG or 
boundary time t coincide, i.e., 

\imv = t. (2.30) 

r->-0 

Now the relevant boundary counterterms are given by [24,31,32]: 

t _ ^divergent . T finite /o oi\ 

l ct — l ct + l ct l^-U 

where 

C er9ent = tttV / d ' x V=7 f 6 + + -V In p (2.32) 

167tG 5 JdM 5 , P =e V 2 12 

+ 1 7^0^-0 + hip), 

llT lte = tt^t f d'x y/^U 4 + S 2 7 ^<W + 5 3 i? 7 2 ) • ( 2 -33) 
1ottG 5 Jqm 5 , p =6 \ J 

These expressions are written using 

7^ dx j dx j = p~ 2 Gij(x k , p) dx { dx j . (2.34) 

Implicitly then, the boundary action contains many potentially divergent factors, e.g., 
y/— 7 = 0(e~ 4 ). Above, we have separated the counterterms which diverge in the limit 
e — > from the finite counterterms. Further W is the Ricci scalar constructed from 
7ij treated as a four-dimensional metric with an external parameter p. Note that even 
though the boundary metric (2.28) is flat, the term proportional to i? 7 2 contributes 
to the holographic renormalization in the presence of space-time varying sources [32]. 
As usual, the coefficients Si of the finite counterterms (2.33) are arbitrary constants, 
which are related to ambiguities in the renormalization scheme. 

Combining all of these expressions to compute the renormalized one-point functions 

for 

(T„>= (e,V,V,V^, (2.35) 
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and (0 3 ), we find 

3 1 1 111 7 

87rG 5 £ = - -a A - -^(p'of ~ ^Poaip' + -p 2 a\ - -p p 2 + ^PoPo + ^ggPo 

_|_ gambiguity 

87rG 5 P = - -a 4 - — (po) + -p aip - — p Oi + qPoP2 - y^oPo + ^Po 

j ^ambiguity 

16vrG 5 (0 3 ) J-pl - Ipg - 2a lP ' + ^ a? - 2p 2 + eg"**** . 



(2.36) 

Here we employ the superscript 1 ambiguity'' to denote renormalization scheme ambi- 
guities introduced by the finite counterterms (2.33) 



pom**** = _ 2W) 2 _ 2§spM _ 1^4 + > (2.37) 

As expected, the physical one-point functions in eq. (2.36) are invariant under the 
residual diffeomorphism (2.22) of the background metric (2.11). 

These one-point correlators are subject to various Ward identities [31]. In particu- 
lar, one has the diffeomorphism Ward identity 

&(T ij ) = -(Gz)d j p . (2.38) 

Of course, when the mass parameter is held constant, this expression reduces to the 
conservation of energy and momentum in the boundary theory. In the present case, 
the configurations are independent of the spatial coordinates and so only the j — t 
component of eq. (2.38) is nontrivial, yielding 

d t £ = (O 3 )d tP0 . (2.39) 

With a time- dependent po, the contribution on the right-hand side describes the work 
done by varying the fermionic mass in the boundary theory. 5 We also note that this 
equation is scheme independent, i.e., eqs. (2.38) and (2.39) hold independent of the 



5 Later we will see that po = V2m,f. 
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values of S^. Finally we observe that the constraint (2.16) reduces to this identity on 
the asymptotic boundary p = 0, i.e., eqs. (2.21) and (2.39) are precisely the same. 
We also have the conformal Ward identify 

(T i i ) = -Po (Os) + — L- QpoDpo + 6 n(j%)) , (2.40) 

where 

S = -S 2 + 65 3 . (2.41) 

Examining the right-hand side of this expression, the first term represents the 'classical' 
contribution from breaking conformal invariance by the introduction of a dimensionful 
coupling p . Note that this term carries an overall factor of (A — d), which happens 
to be -1 in the present case. The remaining terms can be interpreted as anomalous 
contributions. The quadratic term is a standard scalar anomaly that was first studied in 
[33]. In a curved space background, this term would take the form of a Paneitz operator: 
p^d—W /Q)p [34]. The quartic term was previously observed in [24]. The final term is 
a scheme dependent total derivative that is induced by the finite counterterms. While 
the S 2 and 5 3 counterterms in eq. (2.32) are not conformally invariant, they still respect 
scale invariance and so can only contribute with a total derivative in eq. (2.40) [35]. 
Note that with 5 3 = S 2 /6, these two counterterms combine to form a conformally 
invariant combination and hence the total derivative contribution above vanishes. Of 
course, the usual curvature-squared contributions to the trace anomaly [36] do not 
appear on the right-hand side of eq. (2.40) because we are examining the boundary 
theory in flat space, e.g., see eq. (2.28). 

2.1.2 Renormalization of m 2 = —4 bulk scalar 

Constructing the transformation between EF and FG coordinates, using (2.23) for 
m 2 = —4, we find 

v = t-p + 0(p 5 hip), 

1/11 1 \ (2-42) 

r=-[\--p ai (t) + -p 2 a[(t) - -p 3 a[(t) + O(pHnp)) . 

Note that as in the previous case, we have v — > t on the asymptotic boundary, i.e., as 
p ->■ 0. 

In this case, we write the relevant boundary counterterms as [24,31]: 

I ct = Jdivergent + j finite ^ ^ ^ 
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with 



lS° £r9ent = TTtV / d ' x ^( 6 + + oT- <^ > ( 2 - 44 ) 



' dMs,p=€ 



ir te = t^t ^^ r^^ + p^ ■ (2.45) 

lD7rG 5 JdM 5 ,p=e \m p ln P / 

Again we have separated the counterterms which diverge in the limit e — > from those 
which remain finite. Further recall that the coefficients Si of the finite counterterms 
(2.45) are arbitrary constants reflecting the ambiguities in the renormalization scheme. 
Notice that there is a finite counterterm involving the curvature term, i.e., ^ -R 7 0, 
even though the divergent counterterm i? 7 is absent. The absence of such a divergent 
counterterm was established in [24]. Since the bulk action (2.9) is even in 0, all 
divergences must be even in as well. However, the candidate counterterm in question 
is linear in and so does not occur. 

Transforming the expressions in eq. (2.23) to FG form and computing the renor- 
malized one-point functions of (T^) and (0 2 ), we find 

8vrG 5 £ = - ^p l oPo - \a A - ^(ptf - ±pg + S amb ^ , 



where 



16ttG 5 (G 2 ) =~Po + OT U9U * v , 



gambiguity =1^(^)2 ? 



V ambi 9 uity = _ 1^)2 + 5 2 (p l )" , ( 2 ' 47 ) 
pamUguity =2 ^ Q . 

As expected, the physical one-point correlation functions in (2.46) are invariant under 
the residual diffeomorphisms (2.25) of the background metric in eq. (2.11). 

Again, the one-point correlators satisfy Ward identities. In particular, the diffeo- 
morphism Ward identity now becomes 

d l (T tJ ) = -(O 2 )d,p l . (2.48) 

As in the previous case, this equation contains a single nontrivial component here, 6 

d t £ = (0 2 ) d t p l . (2.49) 



'Note that we will see p l = y / 2/3™^ below. 
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As before, this expression is scheme independent, i.e., it holds for arbitrary Si, and it is 
equivalent to the constraint (2.24) when the latter is taken to the asymptotic boundary. 
The conformal Ward identify becomes 

(T t <) = - V (C 2 ) + — Q(^) 2 - 65 2 DtfS) , (2-50) 

which has the same general form as described with the fermionic mass operator in the 
previous subsection. 

2.2 High temperature equilibrium thermodynamics 

Having established the equations of motion for the bulk theory (2.9) of gravity coupled 
to a free scalar, as well as the asymptotic expansions required to extract the one-point 
functions, we apply these results here to examine time-independent or equilibrium 
configurations. High temperature equilibrium thermodynamics of the M = 2* gauge 
theory has been extensively discussed in [23-25]. In the following, we highlight some 
of the salient results of these investigations. In particular, this allows us to connect 
the bulk parameters to those in the M = 2* gauge theory, e.g., as given in eqs. (2.62) 
and (2.63). Note that our nomenclature here should be interpreted as follows: 'high 
temperature' implies that m^j/T <C 1 while 'equilibrium' implies that the bulk gravity 
solution will be time-independent. 

Recall that in the high temperature limit, we are working perturbatively in the 
amplitude of the bulk scalar. Hence the zero'th order solution will be a planar AdS 
black hole with the bulk scalar field set to zero. To leading order in the expansion in 
rribj/T, it is sufficient to solve the linearized equations for the scalar <fi in the black hole 
background. The scalar only backreacts on this background geometry at the quadratic 
order. Thus given the ansatz (2.11), we write 

<j ) (v,r) = X<j ) e 1 (r) + 0(X 3 ), 

A(v, r) = r 2 - ^ + A 2 /i 2 A e 2 (r) + 0(A 4 ) , (2.51) 

H(v,r) = r + X 2 fi Z e 2 (r) + 0(X 4 ) , 

where we have introduced A as an expansion counting parameter. Setting A = leaves 
the AdS black hole, where r = n is the position of the event horizon. Note that we are 
using the superscript e above to denote quantities corresponding to the equilibrium (or 
static) solution. 
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From eq. (2.13), the linearized scalar equation becomes 

r(r 4 -/i 4 ) ° r<Pl + r*-» 
It is convenient to introduce a new radial coordinate p, 7 



= Ol + ^ ^ + Z4— 7 # • ( 2 - 52 ) 



P = - , (2.53) 
r 

such that the asymptotic boundary (r — >■ oo) occurs at p — > while the horizon (r — >• /i) 
is positioned at p — > l~ (to zero'th order in the A expansion). In terms of p, eq. (2.52) 
takes the form 

The solutions are readily expressed in terms of hypergeometric functions and then in 
each case, they take the form: 



For m 2 = —3, 



i? = 7r-^r('?)V2F 1 (|,|,l,l-p 4 )- ( 2 ^) 



This particular solution is chosen by demanding regularity at the horizon. The overall 
normalization constant is chosen here in such a way that 

Po = \», (2.56) 

where po is the leading coefficient in the asymptotic expansion of the scalar given in 
eq. (2.20). Given the above solution (2.55), we can also extract the next independent 
coefficient at order p 3 (oc r~ 3 ) in the latter expansion as 

P2 = -A p 3 -4*- . (2.57) 



Further given the linearized solution (2.55), we can choose the asymptotic solution of 
^4-2 as 

A e 2 (p) = - 1 - + Cp 2 + 0(p 3 \np), (2.58) 
where C is an arbitrary constant. With this choice, we identify in eq. (2.20) 

ai = , a 4 = -p A + \ 2 p 4 C . (2.59) 



7 Note that this radial coordinate is distinct from the the FG radius p, which is related to the EF 
radius r in eqs. (2.29) and (2.42). 
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As noted at eq. (2.22), a x = corresponds to an implicit gauge choice that was made in 
choosing the form of A^p) in eq. (2.58) — compare this expression to the general form 
in eq. (2.20). Similarly, the constant C can be chosen arbitrarily with a gauge choice 
which fixes the position of the horizon, i.e., r 4 = /i 4 (1 — \ 2 C + (9(A 4 )). In a discussion 
of the equilibrium thermodynamics, the most convenient choice is to simply set C = 0. 
However, we will see the analogous constant cannot be avoided in the following studies 
since the quenches naturally cause the horizon to grow. 

Now using eq. (2.36), we can compute the equilibrium expectation values for the 
energy density S e , the pressure V e and the fermionic mass operator {O s ) e : 

8nG 5 S e = (1 - A 2 C) + //A 2 + 0(A 4 ) , 



SvrGs V e = i/i 4 (1 - A 2 C) - ^ 4 A 2 V -SL + 0(A 4 ) , (2-60) 

2 6 7T Z 

16ttG 5 (0 3 > e = /i 3 A 2 J-SL + G(A 3 ) . 

Note that to the order presented, there is no ambiguity in {£ e , V e , ((9 3 ) e }. From 
eq. (2.37), we can see that for a time independent mass, the ambiguities first arise at 
order A 3 in (0 3 } e , while they only arise at order A 4 in £ e and V e . 

Next we compare eq. (2.60) with the energy density and the pressure given in [25] 

8 I 3tt 4 T 2 



4 



, ,x (2-61) 

tt 2 N 2 T 4 1 - 



7T 4 T 2 

for the high-temperature thermodynamics of Af = 2* gauge theory. This comparison 
allows us to identify 

(2.62) 

Combining these expressions with eq. (2.56) yields 

Po = v^m; ^l + O^^ . (2.63) 

Similarly, the equilibrium expectation value of the fermionic mass term in eq. (2.60) 
becomes 

(O 3 r-^0^^m f T'(l + O (^)). (2.64) 
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For m 2 = —4, the linearized solution for the scalar field becomes 



7T 



P 2 2^1 



1 1 1 1 < 

2'2' M ~ P 



(2.65) 



Again, this particular solution is chosen to be regular at the horizon. Here the overall 
normalization constant is chosen in such a way that 



p l = A fi 2 , 



(2.66) 



in the asymptotic expansion of the bulk scalar in eq. (2.23). Further, the second 
independent coefficient in the latter expansion becomes 

Po = -A fi 2 ln(/i/2) , (2.67) 

given the solution (2.65) above. Without loss of generality, we write the asymptotic 
solution of A e 2 as 



A e 2 ( P ) 



C + 



f ^ln2)lnp--ln 2 p)p 2 + 0(p 3 lnp), 



1 2 

54 + 9 

where C is an arbitrary constant. In this case comparing eq. (2.23), we identify 

1 2 , \ , 1 



ai = 0, 



04 



-/1 4 + Ay ( c + 



54 ' 9 ln2lln^--lnV 



(2.68) 



(2.69) 



As discussed above, gauge choices are implicitly at work in setting a± — and fixing a 
final value for C. 

Using eq. (2.46), we can now compute the equilibrium expectation values for the 
energy density £ e , the pressure V e and the bosonic mass operator {O^: 



SvrGs S e = -p 4 + A 2 p 4 



3„ 1, 5 , 5 l l2 
— C + - In // ln2 In 2 2 

2 4 P 18 54 6 



S. 



ambiguity 



8nG 5 V e = V + A 2 /i 4 (--C - - In p + — In 2 
5 2 P P V 2 4 54 
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324 18 



-In 2 2 



^ambiguity ^(A ) j 



16vrG 5 (G 2 ) e = A p 2 In £ + + W - 



(2.70) 



where 



5. 



—A X 2 4 

ambiguity 2 ^ ' 



ambiguity 



1 



:^iA 2 p 4 , 



(2.71) 



& ambiguity — 2 5iA /i 
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Next we compare the energy density and pressure in eq. (2.70) with the analogous 
expressions in [25]: 



£ e = -ir 2 N 2 T 4 ( 1 + ( In — - l) 



A / 9tt 4 T 4 / ' 



where A is an arbitrary scale in the theory. In this way, we identify 



(2.72) 



[2 m 2 h 1 A 

A= V3^' 5l = "2 ln 2 



(2.73) 



Notice that the renormalization scheme ambiguity 8\ is identified with the ambiguity of 
choosing the scale A in eq. (2.72). Combining these expressions with eq. (2.56) yields 



... . / 2 „2 A , „{™>i 



Similarly, the equilibrium expectation value of the fermionic mass term in eq. (2.60) 

becomes _ 

2 N 2 , ttT / „ /m 4 



< ^ = V3^ m ' ln x( 1+0 UJj- < 2 - 75 ' 

3 Holographic mass quenches at high temperatures 

Next we apply our results here to holographic quenches where the mass parameters in 
the dual gauge theory are varied as in eq. (1.3). Let us sketch the general approach 
here for m 2 = — 3: 8 Choosing a specific time-dependent profile for rrif corresponds 
to fixing the asymptotic function po according to eq. (2.63). Given this boundary 
condition, we then numerically solve the scalar field equation (2.13) which allows us 
to determine the sub leading coefficient p2- Note that as in the previous section with 
the high temperature approximation, we are working perturbatively in the amplitude 
of the scalar and so this equation is solved to linear order in A with the static black 
hole background (i.e., setting A = r 2 — fi/r 2 and E = r). Given both p and P2, we can 
evaluate the one-point function (O3) using eq. (2. 36). 9 Further, we can integrate the 



8 This discussion extends to the case m 2 = —4 in an obvious way. 

9 We will make the gauge choice a\ = to simplify this expression, as well as later calculations. 
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boundary constraint (2.21) to determine the metric coefficient a 4 . Having calculated 
the latter, we can evaluate the energy density £ and pressure V (at order A 2 ) using the 
expressions in eq. (2.36). Note then that we are able to evaluate all of the one-point 
functions to leading order in the high temperature approximation without actually 
solving for the time- dependence of bulk metric {i.e., without solving for A 2 {v,r) and 
E 2 (i>, r) below). However, the latter will be required to evaluate non-local probes of 
the quenches [29]. 

As noted above, we are still investigating the time-dependent configurations in the 
high-temperature approximation, rnf^/T <C 1. Hence it is still sufficient to first solve 
the linearized equation for the scalar <fi with the equilibrium black hole background 
and then consider the backreation of the scalar on the geometry at quadratic order. 
We thus organize the perturbative expansion of the ansatz (2.11) with 

<j ) {v,r) = \ ( i )l {v,r) + 0{\ 3 ), 



A{v, r) = r 2 - ?L + A 2 p 2 A 2 {v, r) + 0{\ A ) , (3.1) 



4 



r 

£(u,r) = r + A 2 ii S 2 (v,r) + 0(A 4 ), 

which is a simple extension of eq. (2.51) to a time-dependent situation. Then intro- 
ducing coordinates, 

p = — and r = a v , (3.2) 

r 

the scalar field equation (2.13) becomes, to linear order in A: 

1 3 3 -I- n 4 m 2 

= Ol - o(l - P 4 ) Ol - Z-arfa + d p <h + —,<P l . (3.3) 



Of course, this equation reduces to eq. (2.54) with a time- independent scalar. Further 
at order A 2 , the metric equations (2.14) and (2.15) yield 

=d Tp Z 2 -\{l- P 4 ) 9 2 p S 2 - - 9 T S 2 + - d p ^ 2 + \ d p A 2 + Z 2 --A 2 
2 yy p p 2 p 2 p 

m 2 2 



12p3 ' 

=9%A 2 + I 8 P A 2 + ^ - * 0^ 2 + , 2 - « A 2 



p p- p- p- 



(3.4) 
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while the constraint equations (2.16) and (2.17) become 
=0^2 + 1(1 - P 4 ) 2 - (1 - p 4 ) d%ll 2 - 9 r £ 2 + 9 p £ 2 



+ 5 9tM ~ d ^ d ^ + ^ (9 " 0l)2 + ^ (9r0l)2 ' (3 ' 5) 

2 „ „ 1 



o =^ p e 2 + - «9 p s 2 + — (9^) . 

It is now straightforward to find the asymptotic solutions (as p — > 0) of eqs. (3.3-3.5). 

■ When m 2 = —3 

0i =Pi,o P + p'i, P 2 + (pi,2 + \p% Inp^ p 3 + 0(p 4 lnp) , 

S 2 = - ^p 2 lfi p - ^i )Q pi )0 P 2 + ^(P 3 In p) , (3.6) 

4a = - ^p 2 ,o + ^2,4 + \ ((p'l.o) 2 - Pi,oPi, ) lnp) P 2 + 0(P 3 lnp) . 

In addition, the first equation in eq. (3.5) provides the following constraint: 

1 12 

= 4,4 + 3 (pi,o pi )2 - p'i,o Pi,2) + y^Pi,o Pi,o - gPi,o P'" • ( 3 - 7 ) 

In eqs. (3.6) and (3.7), we have pi i0 = Pi,o( r )) Pi, 2 = Pi,2( r ) and a 2j4 = a 2i 4(r), while 
the prime ' denotes differentiation with respect to r. 
Comparing eqs. (2.20) and (3.6), we identify 

Po = Ap pi i0 , p 2 = Ap 3 (p h2 + p" ^ lnp) , 

( 1 \ (3 ' 8) 

ox = , a 4 = -p 4 + A 2 p 4 ( 02,4 + - ((Pi.o) 2 - Pi,oPi, ) in P J • 

Note that ai = was an implicit gauge choice made in writing £ 2 and A 2 in eq. (3.6). 
Further we note that the boundary constraint (3.7) is precisely the constraint (2.21) 
evaluated at order A 2 . 

As argued above, it is enough to solve numerically only the scalar field equation 
(3.3) to compute p 2 for a specified profile po- Then we obtain a 2j 4 by directly integrating 
the constraint (3.7), 

1 5 2 2 f T 

d2,4 =C 3 ~ -Pl,0 Pl,2 - ^(P'l.o) 2 + gPl,0 Pl,0 + 3 J ds Pl,2(s) Pi, ( S ) , ( 3 - 9 ) 
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(3.10) 



where C3 is an arbitrary integration constant which we will fix to a convenient value 
below. Now we have determined both p 12 and a 2i4 for a specific profile pi >0 (as well as 
having set a± = 0). Thus, as observed above, we have sufficient data to evaluate the 
one-point correlation functions of Ty and O3 using the expressions in eq. (2.36) for a 
mass quench to 0(\ 2 ). 

■ When m? = —4 

0i = - {pifl + p[, In p) p 2 + 0(p 3 In p) , 
£ 2 =0(p 3 \n 2 p), 

M = (a 2A + (^(M,o) 2 " jjpi.oPi.o) Inp - i(M,o) 2 ln 2 p) p 2 + 0(p 3 ln 2 p) • 
The first equation in eq. (3.5) now yields the following boundary constraint: 

= 4,4 " ^Pi,o Pi,o + ^(Po)' Pi,o + ^p!,o (Pi l0 )' + \pi,q Pi,o • (3-H) 

As before, the boundary profiles are functions of r, i.e., p[ = p[ t0 (r), pi t0 = Pi,o(t) 
and a 2j4 = a 2j4 (r) and the prime denotes derivative with respect to r. 
Comparing eqs. (2.23) and (3.10), we identify 

Po = V 2 Pi,o » Po = -Ap 2 (pi,o + p!, ln P) , 

ai = , a 4 = -p 4 + AV ^a 2>4 + ^(pIo) 2 ~ ^Pi,oPi,o^ ln P ~ ^(Pi,o) 2 ln2 P^) ■ 

(3.12) 

Again, note that it suffices to solve numerically only eq. (3.3), which allows us to 
extract p± t o for a given mass profile specified by p\ and then compute the one-point 
function (O2) using eq. (2.46). These expressions also allow us to evaluate (T^-) after 
we determine a 2i4 by integrating eq. (3.11): 

a 2 ,4 = C 2 - ^(pi, ) 2 - \plfl + ^Pi,o p!, ~\j ds (Pi,o(s))' , (3.13) 

where C 2 is another integration constant, for which we make a convenient choice below. 

3.1 Physical observables and ambiguities 

As already discussed, the physical observables which we track during the quenches are 
the one-point correlation functions of the stress-energy tensor and the operator 
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dual to the bulk scalar <fi which induces the quench, i.e., O3 operator for m 2 = —3 
and O2 for m 2 = —4. Further we will consider two types of mass quenches: First, 
we start at t — > —00 with the CFT (i.e., mf^ = 0) in a thermal equilibrium state 
at temperature Tj, and after the quench as t — > +00, it will equilibrate to a thermal 
state of the massive theory with a final temperature Tf but still in the regime where 
rrif h <C Tf. The second scenario will be mass quenches where we start at t — > —00 
with the massive theory in thermal equilibrium with <C Tj and the quench will 
take the mass in the boundary theory to zero. Hence as t — > +00, it equilibrates to a 
new thermal state now of the CFT with temperature Tf. 

3.1.1 m 2 = -3 

From substituting the expressions in eq. (3.8) into eq. (2.36), we find 

3 /3 1 11 1 \ 

87rG 5 E =-/i 4 - //A 2 ( -a 2 , 4 + 1 In// (p' lfi ) 2 + —P?,o + ^PifiPn ~ ^P'l,oPi,o J 

_|_ gambiguity _|_ Q(^\^ 

8nG 5 V - A 2 // Qa 2 ,4 + In/, [~\p%P, fi + ^(K,o) 2 ) + ^(K,o) 2 

- ^i.qPi.2 + ^P'i>i,o) + V ambi9Uity + 0(A 4 ) , 
16ttG 5 (Os) = l -^\ [p% - 4 Pl , 2 - 2 In/, p%) + Q^ 9 ^ + 0(\*) , 



(3.14) 



where: 



gambiguity = K ^y2 ^/ j2 ? 

ambiguity = ^y2 ^ ^2 _ ^ (^2 + ^ ^ > ( 3 .15) 
^ambiguity or 3\ // 

As described above, we consider two classes of quenches: 

■ First, we start with a CFT as t — > — oo (r — > — oo) and quench to a massive theory 
as t — > +00 (r — > +00). This quench is implemented as 

lim p 10 = 0, lim p 10 = l, lim p' 10 = 0. (3.16) 
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Further, the asymptotic values of pi )2 match those in the appropriate thermodynamic 
equilibrium, 

lim p 1)2 = , lim p 1>2 = 4^- . (3.17) 

In particular, the latter is given by comparing to eqs. (2.56) and (2.57). Once the 
profile Pi,o(t) is specified, the full function pi )2 (r) will be determined by numerically 
integrating the scalar field equation (3.3). The metric function a 2 ,4 is then determined 
by eq. (3.9) but we must first fix the integration constant C 3 . It turns out that a 
convenient choice is simply C 3 = 0, which sets a 2 ^(r = —00) = for the present class 
of quenches. Matching the r — > ±00 limits of with the corresponding equilibrium 
results, we identify to leading order 

V2m° f 

H = nTi, A = ^ , m/(r) =m f Pi, (r) . (3.18) 

Further we define two scales, A 2 and A 3 , to be associated with the ambiguities in the 
renormalization scheme, 

5 2 = ^hiA 2 , <5 3 = ^lnA 3 . (3.19) 
Now we can rewrite eq. (3.14) in terms of variables of the boundary gauge theory: 

£ =f ttWZ? (l - ( 2 a 2 , 4 + \{p' lfi f In ^ + \{p' lfi ? - \p^P% 
2 \K) 2 /(m° f ) 4 

V = l -v 2 N 2 T? (l - (2a 2A + l -(p' lfl ) 2 - | Pl>QPl>2 + 2 -p lfi p% (3.20) 

Notice that at the order that we are calculating these quantities, the scheme-dependent 
ambiguities (the arbitrariness in the choice of A 2 and A3) arise only during the evolu- 
tion. That is, asymptotically p[ = = p" and hence the terms involving In A 12 above 
vanish in the initial and final equilibrium configurations. Of course, this is in agree- 
ment with our observation in the previous section that there are no renormalization 
ambiguities in the one-point functions (2.60) at equilibrium. 
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Finally, we compute the final temperature and the final energy density /pressure 
after the quench. Denoting 

r(^ 4 2 f°° 

a ^ = T H+oo a2 < 4(r) = Itf- + 3 J J 8 pM Pi,o( s ) . ( 3 - 21 ) 

and matching the r — > +00 limit of the above expressions (3.20) with the equilibrium 
values in eq. (2.61), we find 

T f _ /r(|) 4 1 \ (m° f ) 2 f(m° f ) 4 \ 

^= i+ (t- 2 ^J^4 +o ( 1 2? l )' (3 - 22) 

It is also interesting to consider the entropy of the gauge theory. For either of the 
asymptotic equilibrium states, we can use the thermodynamic formula: S = (£+V)/T. 
Further, we have £j = 3 Vi for the initial equilibrium of the CFT and so we find 



'"'J 



m° f ) 2 m f{m° f ) 4 



l-o<4^ + ^ ■ (3-23) 



This last expression makes clear that the constant a^° 4 , defined in eq. (3.21), directly 
parameterizes the entropy production of the mass quench. In our simulations, we find 
that af A is always negative — see section 5.1 for details. This result is, of course, in 
agreement with the expectation that the entropy density should always increase. 
Let us consider an adiabatic transition in the mass where 

Kol « 1 (3-24) 
everywhere. We expect (and later confirm numerically) that 

Pi )2 (r) ~ pi )2 (r = +00) p 1;0 (t) = K -f- P1i0 (t) . (3.25) 

That is, for a very slow transition, the system essentially maintains thermodynamic 
equilibrium throughout the process and Pi,2(t) simply tracks pi,o(r). Thus, from 
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eq. (3.21), we have 

(l - 2 J°° ds Pl , (s) j/ lfl ( 8 )j = . (3.26) 



"2,4 



Given this result, eq. (3.23) indicates that the entropy density is constant for slow 
changes of masses. Of course, from eq. (3.22), we see that these adiabatic processes 
where the mass is raised still produce an increase in temperature and energy density 
but a decrease in pressure. This behaviour can be intuitively understood with a quasi- 
particle picture of a system undergoing such an adiabatic transition. 10 In such an 
adiabatic process, the occupation number of energy levels is unchanged but increasing 
mass increases the energy of any given level. Hence the net energy density would 
increase and further, the occupation number of the states after such a quench would 
correspond to a higher effective temperature. 

Notice that since we expect that a^° 4 is always negative, eq. (3.22) indicates that 
the temperature and energy density always increase, even for rapid mass quenches. 
However, the change in the pressure can be either a decrease or an increase. We can 
expect that slow transitions result in the pressure decreasing while rapid quenches, 
with a large entropy production, should result in an increase in the pressure. 

■ For the second type of quenches, we start with a massive theory as t — > —oo 
(r — >■ — oo) and quench to a CFT as t — >■ +oo (r — > +oo). This quench is imple- 
mented as 11 

lim p 10 = 1 , lim p 10 = , lim p' 10 — . (3.27) 
In this case, the asymptotic values of pi t2 become 



r(i) 



4 



lim pi )2 = M^- , lim p h2 = . (3.28) 

7T Z 



The metric function a 2t 4 is again determined by eq. (3.9), however, we have the freedom 
to choose a new value for the integration constant C3. A convenient choice in this case 
turns out to be 

C S = i[pi,oPi l2 ] T= _ 00 = -^-, (3.29) 



10 Of course, a quasi-particle picture is not valid for strongly coupled N — 2* plasma. 
n We use an upper symbol ~ to distinguish quenches of this type. 
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which again ensures that a 2) 4(V = — oo) = 0. Now comparing to eq. (2.61), we match 
the r — > — oo limit of to the corresponding equilibrium results by identifying 

r(|) 4 K) 2 /«) 4 \\ V2m° f 



Svr 4 Tf \ n J ) ' " vrT, ' (3.30) 

m f ( r ) = m ° f p lj0 (r) , 

to leading order. Further we introduce the same scales, A 2 and A 3 , as in eq. (3.19) to 
define the renormalization scheme. Then we rewrite eq. (3.14) in terms of field theory 
parameters: 

/ 3 \ ^ 

£ = 3 -n*NX 4 (l " (2S 2 , 4 + + \{v\ fi f ^ ^ + ^o) 2 - ^ A 

2 \ (m°) 2 /(m ^ 4 



4 



p =^at 2 7? (1 - (2a M + ^p- + i(p;,„) 2 - |p,„p 1>2 + A (3 ' 31) 



0, = - 



Once again, the arbitrariness in the choice of A 2 and A 3 only plays a role during the 
time-dependent portion of the quench — as already pointed out in eq. (2.60), there 
are no scheme dependent ambiguities at equilibrium in quenches of the fermionic mass 
term. 

Finally, we compute the final temperature and the final energy density /pressure 
after the quench. Denoting 

r(^ 4 2 r°° 

a~ 4 = lim a 2 , 4 (r) = + / ds p 1>2 ( S ) p[ fl (s) , (3.32) 

and matching the r — > +00 limit of (3.20) with (2.61) we find 



T,~ 3tt 2 + 2™ 2 ' 1 



~oo 




|=l-(^ + 2^ 4 )^ + 0(^) , (3.33) 

^/_ 1+ /2j_(r_ 

+ 3tt 2 ^° 2 ' 4 
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Turning to the entropy, it is straightforward to show that 



So again the entropy production of these quenches is parameterized by a~ 4 . Further all 
of the general comments made above about the transitions from the CFT to the massive 
theory can be extended to analogous statements about the present class of quenches. 
In particular, a~ 4 vanishes for adiabatic transitions where |p^ | <C 1 everywhere. 

Notice that if Pi,o(t) describes a quantum quench from a CFT to a massive theory, 
then 

Pi,o(t) = 1 -Pi,o(t) (3.35) 

describes a quench from a massive theory to a CFT. Because we are working per- 
turbatively in the amplitude of the scalar field, in this 'reverse' quench, we will also 
find 

r (-) 4 

Pi,2(t) = P2fi(r = +00) - P2,o(t) = f P2,o( T ) ■ ( 3 - 36 ) 

Further, comparing eqs. (3.21) and (3.32), we find 

5£ = a~ 4 . (3.37) 

Hence the entropy production of a quench from the CFT to the massive theory is 
precisely the same as in the reverse quench from the massive theory to the CFT. 
Again, this result occurs because we are working perturbatively in the high temperature 
approximation. 

3.1.2 m 2 = —4 

Substituting eq. (2.46) into eq. (3.12), we find 

87rG 5 £ =^ 4 - fi 4 \ 2 Qa 2i4 - ^Pi,op[, - \(p[,o) 2 ln /^ + |[(pi,o) 2 + \pl,^j 

_|_ gambiguity _|_ Q(^\^ 

8nG 5 V =i/i 4 - /i 4 A 2 Qa 2 , 4 + ^Wi,o + ^(M,o) 2 1* A* " ^(p!o) 2 + ^P?,o) 

_|_ spambiguity _|_ £)^ 4 ^ 

16ttG 5 (0 2 > =/i 2 A (p 1>0 +pi, ln/x) + 0«- 6 « + G(A 3 ) , 

(3.38) 
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where: 

S ambigmty = ^ ^2 ? 

pamWsuity = _ ^ /i 4 A 2 ^2 + ^ ^ ^y/ ? (3.39) 
ambigmty ^ ^ ^ _ 

As before, we consider two classes of quenches: 

■ First, we start with a CFT as t — > — oo (r — >■ — oo) and quench to a massive theory 
as t — > +oo (r — > +oo). This quench is implemented as 

lim p[ = , lim p[ Q = 1 , Urn (p^)' = . (3.40) 

r^ — oo r^+oo r— S-±oo 

In this case, the asymptotic values of pi ;0 are 

lim p 10 = 0, lim pi i0 = -ln2. (3.41) 

r^ — oo t^+oo 

The metric function 02,4 is determined by eq. (3.13) but we must now fix the integration 
constant C2. With some foresight, we choose 



~Pi,o P l i,o + ^(p'i,o) 2 + \p\fi 



= ^i n2 + — + iln 2 2. (3.42) 
,54 81 9 K J 

r=+oo 



Note that with this choice, a 2) 4 does not vanish for either r — > ±00. However, we will 
see that this choice simplifies our description of the entropy production below. Now 
matching the r — > ±00 limits of £ and V with the corresponding equilibrium results, 
we identify to leading order 

P = tT* (l + I , A = \/| ' m ' (T) = K)2 P '' o(T) • (3 ' 43) 

Next we introduce two scales, Ai and A 2 , that are related to the ambiguities in the 
renormalization scheme with 
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Then we can rewrite eq. (3.14) using parameters in the dual field theory: 

3 2Ar2 _ 4 / /2 2^ 10 , 1, , N2l 27rTi 

£ = ^JVT? (l - (-a 2 , 4 - -C 2 - -ri^ - -(p< )> In 

+ >lof + ^,„)^ + («!)), ,,4 5 ) 
/ i \// /2 2^ 26 , 

(Pl,o) - ( 3 fl 2,4 - gC 2 + — P 1>0 Pl,0 



1 2. ,2^/ /8 A 2 (m ^ 2 



) 2 ln^-i^ ) 2 + A P O^ + 0^^ 

Note that in our perturbative expansion in powers of (m°) 2 /T 2 , the pressure is formally 
dominated by the term proportional to A 2 ,. However, this scheme dependent term 
vanishes in an equilibrium configuration where (p\ )" = 0. 

Now we compute the final temperature and the final energy density /pressure after 
the quench. Denoting 



a°° 



2,4 



1 1 f°° 

lim a 2;4 (r) = - - In 2 - - / rfs p 1>0 (s) (pi )0 )'(s) , (3.46) 



and matching the r — > +oo limit of the above results (3.45) with the equilibrium values 
(2.72), we find 



Ti \ 36 6 2,4 J tt 4 T 4 V Tf 



Vi \9 3 Ai 3 2 '7 tt 4 ^ 4 V T i 

Turning to the entropy, since this quench begins with the CFT, we may apply the same 
expression as given in eq. (3.23) to find 



Hence with the judicious choice of the integration constant made in eq. (3.42), we 
find that a^ 4 directly parameterizes the entropy production in these quenches of the 
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bosonic mass. As before, our simulations of the bosonic mass quenches seem to indicate 
that a£° 4 is always negative — see section 5.2 for details. Of course, this matches the 
intuition that the entropy density must always increase. 
If we consider an adiabatic transition where 

l(M,o)'l«l (3-49) 

we expect (and later confirm numerically) that 

Pi,o(r) « pi l0 (r = +00) pi )0 (r) = - In 2 p' 1)0 (r) . (3.50) 

Again, this reflects the expectation that for a very slow transition, the system essentially 
maintains thermodynamic equilibrium throughout the process andpi ) o(r) simply tracks 
p[fi{r). Thus, from eq. (3.21), we have 

a? A ^ - l -\n2 (l-2 jH ' dsp[ fi (s) &,„)'(*))= . (3.51) 

Hence we again find that the entropy density is constant for adiabatic changes of the 
mass. Notice that in this case, from eq. (3.47), these adiabatic processes always pro- 
duce an increase in temperature, however, the sign of change in the energy density and 
the pressure depends on the choice of the renormalization scale A 1; which is needed 
to describe thermodynamics in the massive theory, e.g., see eq. (2.72). Even with a™ 4 
always being negative, similar statements still apply for rapid mass quenches. 



■ The second type of quench, which we consider, starts with the massive theory as 
t — > —00 (r — > —00) and makes a transition to the CFT as t — > +00 (r — > +00). 
These quenches are implemented with 12 

Urn p[ fl = 1 , lim p[ fi = , lim (p[ fi )' = . (3.52) 

T — ^ — OO T^ + OO T— >±00 

In this case, the asymptotic values of pi ;0 become 

lim pifi — — In 2, lim pi t0 = . (3.53) 

t — y — 00 ' r— >+oo ' 

The metric function 02,4 is determined by eq. (3.13) and we must now fix the integration 
constant C 2 . Motivated by the desire for a simple expression for the entropy production, 
we choose 

C 2 = 



1~ ~i ' 

"gPl,0 Pi,o 



= -ln2. (3.54) 

- -=—00 



'Here again, we use ~ to distinguish this second class of quenches. 
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This choice of C 2 again means that 02,4 is nonvanishing for both r — > ±00. Now 
matching the r — > —00 limit of with the corresponding equilibrium results (2.72), 
we identify to leading order 



" = ^ ( 1 - ' A = \H » ^) = « 2 ^lo(r) • (3.55) 



36 tt 4 t; 4 y ' V 3 vr 2 7; 2 ' 

Finally, we introduce two renormalization scales, Ai and A 2 , as in eq. (3.44). Then we 
can rewrite eq. (3.38) in terms of field theory parameters: 

3 9 , r9 _ 4 / (1 1„ 1 1, 10 , 1. , 2nTi 

02 = w < m °> ( p « + ln — + (^J J • 

Finally, we compute the final temperature and the final energy density /pressure 
after the quench. Denoting 

1 1 f°° 

a~ = lim a 2)4 (r) = - In 2 - - / cfe p 1)0 (s) (pi )0 )'(s) , (3.57) 

and matching the r — )■ 00 limit of eq. (3.45) with eq. (2.72), we find 
Tf . / 1 1_\ (m») 4 „ ^(mS 



1 - - + + C 



Tj ~ V36 ' 6 2 'V tt 4 ^ 4 V 7? 



^ V 9 3 A i 3 / * T i 

Given these expressions, it is straightforward to shown that the change in the entropy 
is given by 

Sf -1 h^^+oM^ (3 59) 

Hence we again find that a?? 4 gives a direct measure of the entropy production in these 
quenches. Of course, this simple result (3.59) relies on the judicious choice of C 2 in 
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eq. (3.54). Furthermore all of the general comments made above about the transitions 
from the CFT to the massive theory can be extended to analogous statements about 
the present class of quenches. 

Again, we observe that if p' 10 (r) describes a quench from a CFT to a massive theory, 
then 

p[, (T) = l-p[ t0 (T) (3.60) 

will describe a quench from a massive theory to a CFT. Because our calculations are 
perturbative in the amplitude of the scalar field, we will also have 

Pi,o( T ) = Pi,o( T = +°°) - Pi,o( T ) = - In 2 - p 2 ,o(r) (3.61) 

for this 'reverse' quench. Thus, comparing eqs. (3.46) and (3.57), we again find 

a^ 4 = a£ . (3.62) 

Hence the entropy production of a quench from the CFT to the massive theory is 
precisely the same as in the reverse quench from the massive theory to the CFT. 
Again, this result occurs because we are working perturbatively in the high temperature 
approximation. 

4 Numerical procedure 

As described at the beginning of section 3, the focus of our numerical calculations is 
the linearized scalar wave equation (3.3). The non-normalizable mode of the scalar 
specifies the time-dependent profile of the corresponding mass parameter in the dual 
gauge theory. Solving the scalar equation allows us to extract the normalizable mode, 
with which we can evaluate the one-point function of the corresponding operator. With 
this information, it is also straightforward to integrate the corresponding boundary 
constraint, eq. (3.9) or (3.13), which then allows us to calculate the energy density and 
pressure. 

For the sake of clarity, we describe the numerical procedure for the fermionic and 
bosonic operators, i.e., m 2 = —3 and m 2 = —4, separately. Below and in the following 
section, we also focus on the quenches which go from the CFT to the massive theory. 
As described in the previous section, the reverse quenches going from the massive 
theory to the CFT are then easily obtained using eqs. (3.35) and (3.36) for m 2 = —3 
or eqs. (3.60) and (3.61) for m 2 = —4. 
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4.1 m 2 = -3 



We study mass quenches generated by 



1 1 1 T 1 1 1 t (A l\ 

Pi,o=2 + 2 tanh ~ = 2 + 2 T' ^ ' 

where 

Here we have used the fact that at the boundary t — v — ^ — see eq. (2.30). By 
varying the constant a, we are able to study quenches that are much faster or slower 
than the characteristic thermal time-scale in the corresponding plasma, i.e., l/Tj. 

Solving the linearized equation (3.3) for 1; we are then able extract the normal- 
izable coefficient Pi^ij), which appears as shown in eq. (3.6). In practice, it is more 
convenient to integrate the evolution of 0i(r, p), which we define as 

0i ( T , p)=PPi,o + P 2 drPi,o + ^P 3 In p d* T pi, + p4>i(r,p). (4.3) 

Comparing this definition with eq. (3.6), we see that asymptotically as p — > 0, 

h=Pi,2p 2 + 0(p 3 \np). (4.4) 

Hence organizing the scalar in this way provides for a natural way to enforce the 
desired boundary condition for p 10 and to deal with a regular field in the numerical 
implementation. The scalar wave equation (3.3) now becomes 

o = & T jk - X -^f- d 2 J, - 1 d T ^ + d p ^ + ^ ^ + J , (4.5) 

with 

J = - A p{2 + 3 lnp) c^ rPl ,o + Jjp 4 (2 + 3 lnp) d 2 T p lfi + 2p 3 d T p lfi + X -p 2 p lfi . (4.6) 
From eq. (4.4), the required boundary condition is 0i = 0. 
4.2 m 2 = -4 

We study mass quenches generated by 

Pifi =2 + 2 a = 2 + 2 T ' ( ^ 
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where 

We are again using t — v — - at the boundary and we will vary a to study quenches 
that are much faster or slower than the characteristic thermal time-scale, i.e., 1/Tj. 

Solving the linearized equation (3.3) for 0i, we extract the normalizable coefficient 
Pifi(r), which appears as shown in eq. (3.10). In this case, it is convenient to introduce 
(f>i (r, p) defined as 

0i (r, P) =p 2 lnp p[ fi + p 3 \np d T p[ fi + p fa (r, p) . (4.9) 
Comparing this definition with eq. (3.10), we see that asymptotically 

fa=Pi, p + 0(p). (4.10) 
Writing the scalar equation (3.3) in terms of fa, we have 

o = d 2 p fa - d 2 pp fa - ^ d T fa + i±i^ d p fa - fa + Jo , (4.ii) 

with 



Jo = \p{2 + 3 lnp) d 2 T p[ fi + ^p 4 (2 + 3 lnp) d T p[ fi + 2p 3 (l + lnp) p% . (4.12) 

From eq. (4.10), the boundary condition is 0i = for p — > 0. 

Eqs. (4.5) and (4.11) were discretized employing second order finite difference ap- 
proximations as discussed in Appendix A. The normalizable response, i.e., p\^ for 
m 2 = —3 and pi j0 for m 2 = —4, is then obtained by a straightforward fit from the 
numerical solutions obtained. Examples of the computed values versus time for the 
case a — 1 are shown in fig. 1, together with the adiabatic responses from eqs. (3.25) 
and (3.50) for guidance: 

[P^(r)] admbatic = -^Pi,o(r), (4.13) 
[P^(r)] adiabatlc = -ln2pi,o(r). (4.14) 



5 Results 

In this section we present the results for a range of different quenches with both the 
fermionic and bosonic operators. We discuss each case separately. 
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Figure 1: (Colour online) Evolution of the normalizable component p\^ (left panel) 
and pifl (right panel) during the quenches in eqs. (4.1) and (4.7), respectively, with 
a = 1. The dashed red lines represent the adiabatic response given by eqs. (4.13) and 
(4.14), respectively. 

5.1 Quenches with 3 

The source px,Q is given by eq. (4.1) and hence a indicates the the characteristic time 
scale of the quench. The response of the normalizable component p\^ depends crucially 
on whether the transition is 'fast' with a < 1 or 'slow' with a > 1. Our results for 
both cases are presented in fig. 2. The curves in the left panel present the (rescaled) 
profile of pi t 2 as function of ^ for fast quenches. As indicated, here we have adopted 
values a = {1, 0.5, 0.1, 0.05, 0.01, 0.005}. The right panel we shows the evolution of pi^ 
as a function of - for slow quenches with a = {1, 2.5, 5, 10, 50, 100}. 

Consider first the fast quenches in the left panel of fig. 2. The response of these 
quenches deviate very far from the adiabatic response (4.13), which is represented by 
the red dashed line (with no rescaling). Indeed for the adiabatic response |pi,2| < 0.3 
while the plots show that the maximum of the response p\p scales faster than -\ for 
small a. There are three distinct phases in the evolution of px,i- I n the infinite past the 
system starts at equilibrium with p 10 = and hence P\,i{— oo) = for all quenches. In 
the infinite future, the system ends at equilibrium in the massive theory with p 10 = 1 
and hence pi^+oo) = — r (|) 4 /tt 2 for all quenches. 13 The component of p ii2 which 
scales as 1/a 2 is first excited whenever the source pi i0 noticeably deviates from 0, i.e., 

13 Note that the curves in fig. 2 approach different values at large t because we have rescaled the 
response with a factor of a 2 in this plot. 
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Figure 2: (Colour online) The curves on the left plot represent the evolution of the 
a-rescaled normalizable component, a 2 pi^, as a function of - during the quench (4.1) 
with different values of a. The curves on the right plot represent the evolution of 
the normalizable component, pi^, as a function of ^ during the quench (4.1) for the 
representative values of a. 

whenever p 10 > 0.01 which corresponds to - > —2.3. In fig. 2, a 2 pi^ also appears to 
rapidly approach the asymptotic value 14 for large r when the source is sufficiently close 
to 1, i.e., whenever (1 — Pi,o) < 0.01 and correspondingly - > 2.3. Indeed, most of the 
variation in o?p\^, occurs approximately within the interval - G (—4, 4). In this sense, 
the a 2 -rescaled response follows the characteristic time-scale of the source (4.1). That 
is, faster variations in the source pi$ (smaller values of a) result in both the faster 
excitation and the faster equilibration of a 2 pi,2- 

The right panel in fig. 2 presents the evolution of the normalizable component p±^ 
for slow quenches, i.e., for a > 1. Note that we do not rescale pi^ in this plot. Here 
we see for progressively larger values of a, the profiles come closer to resembling the 
adiabatic response given by eq. (4.13) — the latter is represented by a barely visible 
red dashed line. 

While it is not surprising that in the limit a — > oo the response p\^ follows a 
universal adiabatic profile (4.13), it is remarkable that the response is also universal 
for abrupt quenches, i.e., in the limit a — > 0. In this regard, we first observe that in 
fig. 2, the profile a 2 pi,2 has a relatively simple form for very small a and further 

14 As we discuss below, the curves fig. 2 do not allow us to resolve in the details the approach to the 
final equilibrium. 
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Figure 3: (Colour online) J 7 (a) quantifies the response pi )2 for abrupt quenches, i.e., 
as a — > — see the definition in eq. (5.1). The blue dots correspond to J-{ol^) for 
a n = 2~ n with n = 5, . . . , 11. The dashed red line represents the linear fit to these 
points given in eq. (5.2). 

that this simple form is remarkably similar to a 2 p'{ = — sinh (t/ck)/ cosh 3 (r/a). This 
observation motivated us to define J 7 {a) as the constant which minimizes the following: 

||pi,a - J 7 (a) Pi, || = / d(Pi,o) (pi,2 - ^(«) Pi,o) 2 

/oo 
^Pi.o (Pi,2-^( Q ) Pi,o) 2 • ( 5 - x ) 
-oo 

The first expression defining the norm above was chosen since it seems natural if J-"(a) 
is to describe the universal character of the response. From the second expression, we 
see that since p[ = l/(2a cosh(r/a) 2 , this measure weights most heavily the region in 
the vicinity of r = 0, where the peaks in p 1)2 occur. Fig. 3 shows J-(a n ) for a n = 2~ n 
with n = 5, . . . , 11. As also shown in the figure, these points are very well fit with a 
simple linear expression which takes the form: 

F\ f . t ~ 0.43 - 0.50 In a. (5.2) 

This simple result indicates that the leading contribution in response pi j2 actually grows 
as approximately a~ 2 : ln(l / \fo) as a — > 0, i.e., faster than a~ 2 . Now let us define the 
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Figure 4: (Colour online) Universality of the subtracted response p 12 (defined in 
eq. (5.3)) for abrupt quenches. The different curves are virtually indistinguishable 
from each other. 

'subtracted' response pi j2 as follows 



which is then plotted as a function of in fig. 4. Remarkably, all four curves in this 
plot are almost indistinguishable from each other and so we see that pi^ also has a 
simple universal behaviour in the limit a — > 0. 

Eqs. (3.22) and (3.23) show the characteristics of the final equilibrium state resulting 
from the quench, i.e., the temperature Ty, the energy density Sf, the pressure Vf and 
entropy density Sf, relative to the initial state parameters. All of these ratios depend 
on coefficient a?? A , given in eq. (3.21). As noted in the discussion in section 3.1.1, 
in the adiabatic limit, a — > oo, vanishes and so the entropy density is constant. 
Further for these adiabatic transitions from the thermal state of the CFT plasma to 
the mass-deformed thermal state result in an increase of the temperature Tf > Ti and 
the energy density Sf > Si, while the pressure is decreased Vf < V%. In fact, in all 
of our simulations with finite a, we found that af A < and further jaJfJ grows as 
a becomes smaller — see below. Of course, in agreement with one's intuition then, 
eq. (3.23) indicates that the entropy density increases in a generic quench and in fact, 
the entropy production grows as the quenches become faster. Further, according to 



Pi,2 = a 1 (pi, 2 - J 7 (a) p" ) , 



(5.3) 
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Figure 5: (Colour online) Log-log plot of coefficient ( function of a. The 

dashed red lines represent the linear fits to the data (blue curve) for 'fast' quenches 
with small a (In a — > —00) and for 'slow' quenches with large a (In a — > +00). The 
dashed green and orange horizontal lines indicate the thresholds given in eqs. (5.4) and 
(5.5), respectively. For values of (— af 4 ) above the dashed green line, both classes of 
quenches produce a final energy density which exceeds the initial energy density. For 
values of (— a!f 4 ) above dashed orange line, the final temperature is always larger than 
the initial temperature for either type of quench. See the discussion in the main text. 

eq. (3.22), more rapid quenches result in greater 'heating', i.e., the increase in the 
temperature and the energy density grows as a becomes smaller. 

We will consider these results more quantitatively below but first let us consider 
the properties of quenches which make a transition from a thermal state in the mass- 
deformed gauge theory to a thermal CFT state. As described in section 3.1.1, we can 
make a direct translation of any of the previous results to these 'reverse' quenches 
using eqs. (3.35) and (3.36). For these reverse transitions, the characteristics of the 
final equilibrium relative to the initial state are given by eqs. (3.33) and (3.34) and 
these ratios are controlled by a^° 4 , given in eq. (3.32). However, using the previous 
translation, one finds that, as presented in eq. (3.37): 2°,. In the adiabatic limit, 

these transition from the mass-deformed theory to the CFT state result in 'cooling' of 
the system. That is, a£° 4 vanishes for a — > 00 and while the entropy density is constant, 
the temperature and energy density decrease. However, for very rapid transitions, a|° 4 
becomes large and negative and so eq. (3.33) shows that the quenches are again 'heating' 
the plasma with both Tf > Tj and £f > E{. In fact, we see that the thresholds for 
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increasing the temperature and energy density are slightly different. Specifically, from 
(3.33), £f > Si provided 

-^ 4 >^ 4 f, (5.4) 



while Tf > Tj requires 



2r i^) 4 

-< 4 > -^f- • (5.5) 
As shown in fig. 5, both of these thresholds occur close to a = 1: 

S f >Si a < 1.32 , 

t (5.6) 
7/ > Ti a< 0.86 . 

Fig. 5 presents a log-log plot of the coefficient function of a. The behaviour 

of the plot reveals simple behaviours in both the 'slow' and 'fast' regimes with a ^> 1 
and a < 1, respectively. In fig. 5, the dashed red lines show the linear fits to the data 
for these ranges of a. In particular, we find 

slow: ]n(-aS° 4 )L ~ -2.465 - 1.02 In a, a > 1 , 

' lftt (5.7) 
fast: \a(-af A )\ f . t ~ -2.170- 2.02 In a, a < 1 . 

Hence we find that a^° 4 vanishes as a — > oo with a^° 4 oc 1/a. However, we also find 
that a^° 4 is divergent for a — > with a~ 4 oc 1/a 2 . Let us further observe that the 



constant term in the 'fast' fit seems to match: In 
our numerical results seems to indicate that the 
by 



r(!) 4 /(2- 



2^ 



-2.169- ••. Hence 



eading behaviour as a — > is given 



rm 4 1 

which, of course, calls for an analytic derivation. Unfortunately, at present, we can 
provide no insight into such a derivation. However, we note that the scaling for slow 
transitions can, in principle, be deduced analytically from the evolution of <pi with 
a»l — see the discussion below at eq. (6.20). 

Using eq. (3.22), we can translate the asymptotic behaviour in eq. (5.7) into 

AS _£ f - £f iabatic { I for a > 1 , 

t% tl [ ^ -^r- for a < 1 , 

where gjdiabatic j g ^ ne final energy density for an adiabatic transition, i.e., 



Sf abatic = £i ( 1 + + ( ^r" ) ) • (5-10) 
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Note that the same expression (5.9) also applies for AS /Si in the 'reverse' quenches 
describing transitions from the mass-deformed gauge theory to the CFT. Further we can 
similarly deduce the asymptotic behaviour for the relative change in the temperature 
and the pressure using eq. (3.22) and the entropy using eq. (3.23). Using eq. (5.8), we 
have a more precise expression for the limit a — > 0, 

*£..r(j)' ("»?)' i 

£, ~ tt" Tf a* + ' (5 - U > 
In particular then, we note that a discontinuous quench with a = would seem to 
produce physical divergences — see further discussion in section 6. 

However, let us observe that previously we found that the leading scaling of am- 
plitude of the response was a~ 2 ln(l/o;) — see discussion around eq. (5.2). Hence one 
might have expected an even stronger divergence in eq. (5.11), which 'only' scales like 
a~ 2 . This behaviour can be understood as follows: The divergence in eq. (5.11) simply 
reflects the divergence in a^ 4 given in eq. (5.8). Now in eq. (3.21) which defines the 
latter, the integral dominates for small a and so from eq. (5.3), we have 



J ds (V(a) p'{ fi {s) - pi, 2 (s)^ p' lfi 

dspi, 2 (s)p' 10 (s). (5.12) 



3 a 5 

That is, the specific form of the leading response, i.e., being proportional to p'( , leads 
to a total derivative which yields a vanishing contribution above. Hence, only the 
'subtracted' response (5.3) contributes in the integral to produce the a~ 2 scaling in 
eq. (5.11). 

We conclude this section with more quantitative analysis of the excitation and the 
equilibration time scales in the evolution of the normalizable component Pi,2(t). As 
described above, for an adiabatic transition, the system essentially maintains ther- 
modynamic equilibrium throughout the process and Pi^ij) simply tracks Pifiij), as 
indicated in eq. (3.25) or (4.13). Hence, we define a measure of the deviation from the 
local equilibrium as follows 

Pi )2 (r) - bi, 2 (r)] a&: ,„ ; , 



equilibrium 
Pi, 2 



(5.13) 



where [pi, 2 (r)] , , tic is given in eq. (4.13) and 



r 

equilibrium^ ^ p t T \ = . ( 5 .1 4 ) 



t— >+oc 7T Z 
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Figure 6: (Colour online) Extraction of the excitation/equilibration rates for a = 1 
quench. The horizontal green line is the threshold for excitation/equilibration which 
we define to be 5% away from local equilibrium as determined by 8 neq , as given defined 
in eq. (5.13). The dashed red lines indicate the earliest and latest times of crossing 
this threshold, which we denote as r ex (for excitation time) and r eq (for equilibration 
time), respectively. 

Fig. 6 shows a typical plot of 5 neq . We further define the excitation time r ex as the 
earliest time in the quench at which 5 neq when the latter quantity exceeds a predefined 
threshold value, which we choose to be e = 0.05, i.e., 

Sneq < e = 0.05, for all t < T ex . (5.15) 

While the choice of the threshold here is somewhat arbitrary, we do not expect that 
it will greatly effect the following observations. Similarly, we define the equilibration 
time r eq as the latest time at which S neq exceeds the same threshold value, i.e., 

$neq < e = 0.05 , for all t > r eq . (5.16) 

All of these quantities are identified in the sample plot shown in fig. 6. In fig. 7, we 
show r ex and r eq for quenches across a wide range of a. 

Examining fig. 6, one sees that S neq has several the zeros (at finite r/ot) which re- 
flect the fact that pi^ is oscillating about [pi^adiabatic — later we will see that these 
oscillations are described by quasinormal modes of the scalar. The amplitude of os- 
cillations become smaller as a grows and in fact, for a sufficiently large, 5 neq will be 
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Figure 7: (Colour online) Excitation time r ex (left panel) and equilibration time r eg 
(right panel) as a function of a. The threshold is fixed to be e = 0.05 — see eqs. (5.15) 
and (5.16). In both simple linear fit describes the behaviour for small a, as 

shown with the red dashed lines. 

below any fixed threshold e throughout the quench. While we did not try to pinpoint 
the precise value of a at which this occurs, we did establish that this behaviour sets in 
before a = 16 with e = .05. Further, we see in fig. 7 that r ex and r eq behave irregularly 
for large a. In particular, r ex is clearly discontinuous for a ~ 1.8 and there is a similar 
discontinuity in r eq , albeit less profound. Hence, as defined above, these time scales 
are not useful diagnostics for slow quenches. 

On the other hand, both r ex and r eq exhibit simple behaviours for small a, as shown 
with the linear fits in fig. 7. In the left panel, the linear fit for the excitation time is 

~ 2.8-1.2 ma. (5.17) 

fit 



a 



A heuristic understanding of this behaviour is as follows: In the discussion around 
eqs. (5.1-5.3), we demonstrated that the leading response takes the form 

Pi.2(r ex ) ~ F{a)p%{r ex ) ~ e 2 ^ Z« . (5 .i 8 ) 

cr a 

where the latter uses the approximate form of p'( Q for large negative r. Hence if we set 
a fixed threshold for the response at early times, we should expect that 

Pi2(r ex ) ~ const lna + C(ln(-ln(a))) , (5.19) 

a 

which matches the result given in eq. (5.17). Using eqs. (4.1) and (4.2), this behaviour 
can be translated to an excitation time in terms of the original boundary time 

t ex = T— - T In (l/T/T) for TfT < 1 . (5.20) 

a 
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The scaling of the equilibration time at small a is different. For the right panel on 
fig. 7, the linear fit was 

~ f.f - f.O In a. (5.21) 



ln^ 

a 



fit 



Hence we have for the equilibration time in fast quenches 



^ ~ -, (5.22) 

a a 

or, in terms of the original boundary time, 

t eq = T^ ~ 1 (5.23) 
a Li 

That is, irrespective of how quickly the system is driven away from the initial equi- 
librium, the return to the final equilibrium is determined universally by the typical 
thermal time-scale. 

This behaviour seems closely connected to the fact that the late-time response 
is controlled by the quasinormal modes of the scalar field [37,38]. The quasinormal 
oscillations are well illustrated in fig. 8, which shows the late-time evolution of pi^ 
for a range of small values of a. For all of these fast quenches, a single mode clearly 
determines the decay of pi ;2 after a very short time. A rough fit of the oscillation periods 
and the slopes indicates that in all cases, the decay is governed by a quasinormal 
mode with u/{2'KT i ) ~ (1.095 + i0.87). This value is consistent with the expected 
frequency (1.099 + % 0.879), found in [38]. Beyond this fundamental mode, the first 
overtone can also be extracted for about one oscillation period and a rough fit yields 
real and imaginary parts of its frequency which agree with the expected values to 
within approximately 25%. 

5.2 Quenches with 2 

Our discussion of the results for quenching the bosonic mass operator will be shorter, 
as it parallels the previous discussion of quenches of the fermionic mass operator. 
However, we will see that the physical characteristic of these types of quenches differ 
in many respects. For the present quenches, the profile of the source p\ is given by 
eq. (4.7). Again, the characteristics of response in the normalizable component p^o 
depend on whether the transition is 'fast' with a < 1 or 'slow' with a > 1. Our results 
for both cases are presented in fig. 9. The curves in the left and right panels present pi^ 
as function of - for fast and slow quenches, respectively. We present fast quenches with 
a = {1, 0.5, 0.1, 0.05, 0.01, 0.005} and slow quenches with a = {1, 2.5, 5, 10, 50, 100}. 
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Figure 8: (Colour online) Behaviour of the response coefficient versus time for different 
fast quenches. Note that on the vertical axis, we have defined [pi )2 ] = pi i 2 + r(3/4) 4 7r -2 . 
In all cases, the same quasinormal mode governs the dynamics very soon after the 
quench. 

First we consider the slow quenches, i.e., with a > 1 in the right panel. As found 
before with quenches of O3 operator, progressively larger values of a result in profiles 
that more closely resemble the adiabatic response given here by eq. (4.14). The latter 
profile 'appears' in the plot as a black dotted line, but for the most part, it is covered 
by the profiles for a = 50 and 100 which track this adiabatic response. 

Now we turn to the fast quenches in the left panel. The response shown there 
is very far from the adiabatic curve, which corresponds to eq. (4.14). Although the 
maximum amplitude of the response grows as a becomes smaller and the quenches 
become faster, comparing to fig. 2, we see that the present increase is not as nearly 
as dramatic as found for quenches of O3 in the previous section. Notice that neither 
the response pi t o nor the time r is scaled for fast quenches in fig. 9, in contrast to the 
corresponding plot in fig. 2. With this slower grow, there is no distinctive profile that 
emerges here for p^o in this limit. Hence to determine if the response is scaling with 
a, we first simply consider the peak value of this normalizable component, i.e., 

To(a) = max \pi,o\ ■ (5-24) 

t£( — 00, +00) 
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Figure 9: (Colour online) Normalizable component pi.o as a function of r for fast and 
slow quenches is shown in the left and the right panels, respectively. 

The left panel in fig. 10 shows J-o(«n) for the a select set of values, a n = {1/32, 1/64, 
1/128, 0.0025, 0.0001, 0.00001}. As the red dashed line in the figure shows, these points 
are quite well fit with a simple linear expression which takes the form 

JoL,- -0.80 - 0.94 In a. (5.25) 

The rescaled response, p 10 = — p 1<0 / ^(a), is then shown in the left panel of fig. 11 
as a function of p[ . It seems that these curves are slowly converging to a straight 
line such that in the limit a — > 0, the profile of this rescaled response is simply given 
by Pi,o = Pi 0' which is represented by a dotted orange line in the figure. That is, in 
the limit of very fast quenches, the rescaled response simply tracks the source. We 
must clarify this claim, which seems to suggest that profile of p^o should approach a 
simple step function as a — > 0. Clearly, the latter is at odds with the profiles shown in 
fig. 9. However, note that in fig. 11, we are plotting rescaled response as a function of 
p[ o and with the profile in eq. (4.7), this source varies across 90% of its range within 
— 3a/2 < t < 3a/2. Hence in this plot, we are really focussing the behaviour of pi 5 o 
in a very narrow range around r = when a is small. 15 Further, this observation 

15 To see more clearly that as a — > 0, p± t o is approaching this step function in time, one can plot 
pi,o as a function of r/a - rather than just tau, as was done in fig. 9. However, as is evident from 
fig. 11, the response still falls off rapidly for r/a > 3/2 for the values of a shown there. 
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Figure 10: (Colour online) J- (a) quantifies the response pi i0 for abrupt quenches, i.e., 
as a — >■ 0. The red dots correspond to 7o( a n) (defined in eq. (5.24)) for a n = {1/32, 
1/64, 1/128, 0.0025, 0.0001, 0.00001}. The dashed red line represents the linear fit 
to these points given in eq. (5.25). The blue dots correspond to J-\{a. n ) (defined in 
eq. (5.26)) for the same a n . The dashed blue line represents the linear fit to the latter 
points given in eq. (5.27). 



motivated us to define T\(ot) as the constant which minimizes the following norm: 

||pi,o - F\(a) Pi.oll = / d (Pi,o) (Pi,o - F\(a) p[, Q ) 2 

Jo 




dr(p[, y (Pi,o-^i(a)pi, ) 2 > (5-26) 



in analogy with eq. (5.1) for the dimension-3 operator. As before, the factor (p[ )' = 
(2a cosh(r/a;) 2 )~ 1 in the last expression weights most heavily the integration in a nar- 
row range around r = for small a. The right panel in fig. 10 shows J-~i(oi n ) for the 
the same values a n used above. The blue dashed line in the figure indicates a simple 
linear fit to these points, which takes the form 

Tx | fit ^ -0.54- 1.03 In a. (5.27) 

In this case, the rescaled response, pi^ = — pi i o/^ r i(tt), is then shown in the right panel 
of fig. 11 and again the results suggest a slow convergence towards pi 7 o — p[ Q as a — >■ 0. 
Hence, with both definitions of the scaling factor in eqs. (5.25) and (5.26), the results 
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Figure 11: (Colour online) Limiting behaviour of the rescaled response for abrupt 
quenches with pi t0 = — pi t o/J r o(a) in the left panel and p^o = — pifi/J-i(a) in the right 
panel. In both cases, the curves suggest a slow convergence towards straight line where 
the rescaled response simply follows the source, e.g., pi )0 = p[ on the left, as a — > 0. 

suggest that the same scaling of the 'prompt' response with pi tQ ~ \ogap[ for fast 
quenches, i.e., with small a. 

The reduced response of the present quenches with 2 will also become evident 
below when we consider the properties of the final equilibrium state. Finally, we re- 
iterate that neither pip nor r is scaled in fig. 9 for the fast quenches. Hence we observe 
that the relaxation of the system after a fast quench is almost 'universal', i.e., it is 
essentially independent of a. 

We now turn to eqs. (3.47) and (3.48), which give the characteristics of the final 
equilibrium state relative to those of the initial equilibrium, i.e., Tf/T i7 £f/£i, Vf/Vi 
and Sf/Si. All of these ratios depend on coefficient a£° 4 given in eq. (3.46). In par- 
ticular, this coefficient was defined to give a direct measure of the entropy production 
during the quench, as shown in eq. (3.48). This coefficient is shown for a wide range 
of a in fig. 12. One simple observation is that in our numerical simulations, we always 
found a^ 4 < 0, which must intuitively be the case from eq. (3.48) in order that Sf > Si. 
We might also observe that |a^° 4 | is becoming very small for large a, which agrees with 
the discussion in section 3.1 where we argued that af 4 — > for a — > oo. Let us add here 
that according to eq. (3.47), we may describe the quenches going from the CFT to the 
mass-deformed gauge theory as always 'heating' the system. That is, we always have 
Tf > Ti, even in the adiabatic system. On the other hand, according to eq. (3.58) for 
'reverse' quenches going from the mass-deformed gauge theory to the CFT, the system 
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Figure 12: (Colour online) Coefficient function of a. The dashed red lines 

represent the linear fits to the data (blue points) for fast (left panel) and slow (right 
panel) quenches. Values of (— a^° 4 ) above dashed orange line in the right panel always 
result in quenches where the final temperature exceeds the initial temperature. 

'cools' if the transition is sufficiently slow, sufficiently small. However, as 

the quenches become faster, a^° 4 becomes large and negative resulting in Tf > T again. 
Specifically, from (3.58), Tf > T requires 16 

-a- > 1 . (5.28) 

As shown in the right panel of fig. 12, this threshold occurs close to a = 1: 

T f > T a < 0.58 . (5.29) 

In certain respects, the results for a^° 4 in fig. 12 for quenches induced by the bosonic 
mass term are similar to those in fig. 5 for quenches induced by the fermionic mass 
term. In particular, both figures show a simple behaviour for both a ^> 1 and a C 1 
and a clear transition between these two around a ~ 1. However, there is one striking 
difference in that the fast quenches are shown in the left panel of fig. 12 which shows 
— af A (rather than ln(— a^° 4 )) as a function of In a. The simple fits for the fast and slow 
quenches then take the from 

slow: ln(-aS° 4 )L ~ -2.688 - 1.00 In a, a>l, 

' h (5.30) 

fast: -af 4 L.~ 0.05 - 0.17 In a, a < 1 . 

i i j%t 

16 Here we are using eq. (3.62) which states that a£° 4 = a£° 4 for the reverse quenches described by 
eqs. (3.60) and (3.61). 
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Using eq. (3.47), we can translate the asymptotic behaviour in eq. (5.30) into 

AS £ f - Sf iabatic J \ ^ for a > 1 , 

^ \ 1 t m o^4 (5.31) 



Si Si i wi i n \ (thI 



ln(l/a) ^ fora«l, 



where 

Sr abaUc = Si ^1 - - In + O ) ) . (5.32) 

Of course, we can also determine the asymptotic behaviour for the relative change in 
the temperature and the pressure using eq. (3.47) and the entropy using eq. (3.48). In 
particular then, we note that a discontinuous quench with a = would again seem to 
produce physical divergences. However, we are again seeing a weaker response here for 
the quenches induced by 2 - Recall that the analogous divergence in eq. (5.11) was 
proportional to 1/a 2 for the quenches with (9 3 . 

Further, we note that the scaling of a^° 4 for a — > shown in eq. (5.30) is slower 
than the scaling of the amplitude of the response for the fast quenches found in our 
previous discussion — this was also found for the analogous results with O3. However, 
if we assume the 'prompt' response takes the form pi ; o — logap^ as discussed above, 
then the arguments given for the O3 case around eq. (5.12) would in fact indicate that 
both a^° 4 and p±fi scale in the same way. Hence the reasons for the discrepancy in the 
scaling of these two quantities must be more subtle here than before. 

We now turn to the excitation and the equilibration time scales in the evolution of 
the normalizable component p 10 (r). By analogy to eq. (5.13), we define 



Pi,o( T ) ~ \Pi,o( T )] 



adiabatic 



equilibrium 
Pl,0 



(5.33) 



where [pi >0 (r)] adiabatic is given in eq. (4.14) and 

equilibrium^ ^ ^ ( r ) = _ ln 2 . ( 5 .34) 

t^+oo 

A typical behaviour of S neq is presented in fig. 13. We define the excitation and equi- 
libration times r ex and r eq precisely as before in eqs. (5.15) and (5.16). We indicate 
these times in fig. 13, as well as the threshold e = 0.05. We extracted r ex and r eq for 
a wide range of a and our results are shown in fig. 14. As in the previous discussion, 
we only expect these time scales to be a useful diagnostic of fast quenches with a < 1. 
For example, while we did not pinpoint the precise value of a such that 5 neq < e for all 
r G (—00, +00), we did establish that such a value of a occurs prior to a — 16. 
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Figure 13: (Colour online) Extraction of the excitation/equilibration rates for a = 1 
quench. The horizontal green line is the threshold for excitation/equilibration which 
we define to be 5% away from local equilibrium as determined by 5 neq , as given defined 
in eq. (5.33). The dashed red lines indicate the earliest and latest times of crossing 
this threshold, which we denote as r ex (for excitation time) and r eq (for equilibration 
time), respectively. 

The left panel of fig. 14 suggests that, in the limit In a — > — oo, 

— ~ln(-lna), (5.35) 
a 

however, we must say that our results in this regime were not accurate enough to 
definitively fix this scaling. According to eq. (4.7), the scaling (5.35) would translate 
to 

t ex ~ T lnln(l/T;T) for a < 1 . (5.36) 
In the right panel of fig. 14, we see the scaling of the equilibration time for small a is 



ln^ 

a 



~ 1.5 - 1.0 lna. (5.37) 

fit 



Hence for very fast quenches, r eq /a oc 1/a or in terms of the original boundary time, 
we have 

t eq = T T -^ ex i (5.38) 
a Li 

Much like for quenches of O3, here we find that irrespective of how quickly the system 
was driven away from the initial equilibrium, the return to the final equilibrium is 
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Figure 14: (Colour online) Excitation time r ex (left panel) and equilibration time r eq 
(right panel) as a function of a. The threshold is fixed to be e = 0.05 — see eqs. (5.15) 
and (5.16). In the left panel, it seems like \r ex \/a ~ ln(— In a) in the limit In a; — > — oo. 
In the right panel, the behaviour of ln(r eg /a) is linear in In a for small a, as shown 
with the red dashed line. 

determined universally by the thermal time scale 1/Tj. As with the previous quenches, 
here again the response is essentially controlled by a single quasinormal mode of the 
scalar field at a short time after the quench. This behaviour is well illustrated in fig. 15, 
which displays the evolution of \pi,o] = pifi—p e i^ lhbrium = px t o+la(2). A rough fit to the 
oscillation periods and slopes indicates that the corresponding quasinormal mode has 
u/(27rTi) ~ (0.64+z 0.4), which is consistent with the expected value of (0.644-H 0.411) 
found in [38] . As in the previous section, the first overtone can be extracted for about 
one oscillation period and once again, a rough fit yields real and imaginary parts of 
the quasinormal frequency which agree to within 25% of the expected ones. 



6 Conclusions 

In this paper, we initiated a program to study quantum quenches in strongly coupled 
quantum field theories using holography. We used gauge-gravity correspondence [6] 
to translate the problem of quenching the coupling Aa of a relevant operator Oa in 
planar M = 4 SYM, as given in eq. (1.2), to a classical problem in Einstein gravity 
coupled to a massive scalar field in asymptotically AdS^ spacetime. Here, we focused 
on two cases where A = 2 and 3 but our analysis is readily extended to considering 
general values of A [29]. Our discussion also only considered 'thermal quenches' where 
the system begins in a thermal state of the gauge theory plasma and was limited to a 
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Figure 15: (Colour online) Behaviour of the response coefficients versus time for rep- 
resentative fast quenches. As is evident in the picture, the same quasinormal mode 
governs the dynamics very quickly after the quench. 

high temperature limit. That is, during the quench, A a is always small compared to 
the temperature of the initial state Tj, i.e., 



The latter means that our calculations were perturbative in the coupling Aa- However, 
there was no restriction on the time scale T governing the rate of change of the coupling 
in the transition. In particular, we considered profiles of the form given in eq. (1.3) 
and our results were described in terms of the dimensionless parameter: a = 7rTj T. 
In our analysis, we considered arbitrary values of a and in particular, we examined 
the limits of adiabatic transitions with a — > oo and abrupt quenches with a — > 0. As 
indicated in eq. (1.3), we discussed both the quenches from the conformal SYM theory 
to the mass-deformed gauge theory (with the plus sign for which Aa(£ — > — oo) — > 0) 
and the reverse quenches from the mass-deformed gauge theory to the conformal gauge 
theory (with the minus sign). In all of these quenches, we computed, to leading order 
in the coupling, the response of some basic gauge invariant observables, namely, the 
one-point correlators of the stress tensor (T^-) and the mass operator (0a)- A detailed 
discussion of the results is given in section 5. 



A A | < % 



(6.1) 
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One of the most interesting results coming from our analysis is that the 'infinitely 
fast' quenches seem to be ill-defined because of divergences appearing in the limit 
a — > 0. These divergences appear most dramatic for quenches of O3, the fermionic 
mass operator, in expressions such as eqs. (5.9) or (5.11) which indicate that the energy 
of the final state diverges as I /a 2 as a — > 0. This result is of particular note since 
it is precisely such 'infinitely fast' quenches are studied in the seminal work on the 
topic [4,5]. In their description of these quenches, the evolution from t = —00 to 
0" is essentially regarded as preparing the system in the ground state of the initial 
Hamiltonian and then this ground state is used as the initial condition at t = + 
as the system evolves forward with the new 'quenched' Hamiltonian. In our more 
physical approach to describing fast quenches, the transition from the initial to the 
final Hamiltonian is accompanied by an enormous response in the infinitesimal interval 
around t — 0, as can be seen in the scaling of ^2,0 for O3 and of pi )0 for O2 as a — > 0. 
Hence one might worry that our results call into question the approach of these early 
studies. Of course, an important difference between these studies and our work is that 
the former considered free or weakly coupled field theories whereas we are studying a 
strongly coupled field theory. 

However, we must insert a word of caution with respect to our discussion above. 
The framework considered here treats the bulk scalar field perturbatively. In particular, 
we scaled the amplitude of the field by a factor A, which was then gave a perturbative 
expansion of the results for the boundary theory in terms of m°/Tj or (m") 2 /Tf, using 
eq. (3.18) or (3.43), respectively. However, this approach only properly accounts for 
the backreaction of the bulk scalar on the background spacetime as long as the quench 
is not too fast. Looking at the gravitational equations, we find terms containing time- 
derivatives of the scalar, i.e., terms proportional to (d v <f)) 2 in eq. (2.16), as well as to 
d v (j)d r (j) in both eqs. (2.15) and (2.16). Hence heuristically we should say that in order 
to ensure the full backreaction of the bulk scalar remains small, we must in fact fix 
\/a <^1. In other words, our perturbative calculations only yield reliable results, e.g., 

in eqs. (5.9) or (5.11), as long as f-^f J remains small. This issue is also evident 
in the observation that with a very high final temperature (in particular, much larger 
than the mass scale of the operator in the quench), Tf/T i: Sf/Si and Sf/Si can not all 
diverge in precisely the same way, as would naively be the case with our perturbative 
treatment. Rather, the divergences (or scalings) must arrange themselves such that 
Sf/Si oc (Tf/Ti) 4 and Sf/Si °c (Tf/Ti) 3 . Hence it is clear that our present analysis does 



56 



not provide the final story and to properly understand the behaviour of the 'infinitely 
fast' quenches, we must turn to a full nonlinear analysis of the bulk gravitational 
equations of motion (2.13-2.17) [29]. However, to close this discussion, let us note that 
the nonlinear analysis of [13] do yield divergent results in the limit a — > 0. 

Further, we wish to stress that the above discussion does not invalidate the various 
universal scaling properties and profiles that were found in section 5 for a — > 0. Rather 
these should be interpreted in terms of choosing A <C a m i n for some small value a m i n . 
Then for such a sufficiently small amplitude, the universal behaviours will still correctly 
describe the physics of quenches in the regime a min < a < 1. 

Allowing the couplings of the boundary gauge theory to vary in space and time 
introduces a variety of new ultraviolet divergences into the quantum field theory. The 
power of holography is that it provides a well-defined framework for the renormaliza- 
tion of the boundary theory, even in the presence of these space-time varying couplings. 
As discussed in section 2.1, the bulk gravity theory provides a natural geometric con- 
struction which regulates the new divergences and allows us to identify the appropriate 
counterterms to renormalize them in a covariant way. Recall that roughly, the expec- 
tation value of the quenching operator Oa and its coupling Aa are encoded as normal- 
izable and non-normalizable coefficients (correspondingly) of the dual bulk scalar, 17 
0a = 0a(£,p), 

0a ~ A A p 4 ~ A + O aP a for 2 < A < 4 , 

(6.2) 

02 ~ A 2 p 2 In p + 2 p 2 for A = 2 , 

where p is the usual Fefferman-Graham coordinate (2.26) in the asymptotically AdS 
spacetime. In particular, examining eq. (2.32) for the fermionic mass operator, we 
see that there are additional counterterms proportional to (dX^) 2 and i2A§, which 
play a role in backgrounds where A3 varies in time (and/or space). We note that 
this is completely analogous to placing the boundary theory in a curved background 
spacetime, which creates the necessity of including additional counterterms which are 
functionals of the background curvatures [40] . We also note that some of the new UV 
divergences, e.g., those associated with the two counterterms above, are logarithmic in 
the cut-off scale and as a result, new scheme dependent ambiguities are inevitable in 
the renormalization procedure — we return to this issue below. 

17 We focus here on relevant operators (with A < 4). We can also consider operators with 1 < A < 2 
using the 'alternate quantization' of the holographic theory [39]. For these, the powers of p associated 
with the two coefficients are interchanged in eq. (6.2). 
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Through Einstein's equations (2.7), the bulk scalar backreacts on the geometry at 
least at the quadratic order and hence the renormalized boundary stress-energy tensor 
must depend on the (non-)normalizable coefficients of <p A at least quadratically. Simple 
dimensional arguments imply that for a constant coupling the boundary stress-energy 
tensor (having dimension four) can always depend on the bilinear AaC^a f° r any A. 
This combination is seen to explicitly appear in the energy tensor and the pressure 
given in eqs. (2.36) and (2.46) for for A = 3 and 2, respectively. Continuing these 
dimensional arguments for constant coupling, the stress tensor can also depend on 

4 

A A n for A = 4 - - with n = 2, 3, 4, • • • , 

4 n (6-3) 

Or for A = — with m = 2, 3, 4 . 

m 

For A = 2 and 3, the dependence as in (6.3) was indeed established in [24,25] and can 
again be explicitly seen in eqs. (2.36) and (2.46). 

Once Aa depends on time, additional combinations involving time- derivatives of Aa 
or O a can appear in boundary stress tensor, as long as the dependence has the correct 
mass- dimension four (and is at least quadratic in these coefficients). In particular, the 
stress tensor may depend on 

Aa " 2 d t \ A d t \ A , A A n_1 d 2 tt \ A for A = 4 - - with n = 2, 3, 4, • • • , 

n (6-4) 

d t O A d t O A , O a d 2 tt O A for A = 1 . 

Here we have restricted our attention on relevant operators (with A < 4). We also 
only consider the possibility of an even number of time-derivatives because otherwise 
the corresponding operator could not be extended to a covariant expression. Again, 
these additional terms are explicitly seen in eq. (2.36) for the A = 3 operator. Further, 
as indicated by eq. (6.4), no such contributions appear in eq. (2.46) for A = 2. 

Returning to the possibility that rapid quenches may produce singular behaviour, 
the additional contributions to the boundary stress-energy tensor involving time deriva- 
tives of the couplings Aa, as in eq. (6.4), should be worrisome. In particular, these 
terms suggest that instantaneous changes in the corresponding couplings will produce 
a singular response in the stress tensor. Alternatively, if such a singularity is to be 
avoided, the response of the normalizable coefficient O a should be 'finely correlated' 
to cancel the potential divergence, e.g., through the contributions of the AaC^a terms. 
While in section 5, we saw that for both operators the specific profile of the leading 
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response was closely related to that for the coupling, e.g., see eq. (5.1), we must re- 
iterate that our perturbative calculations can not fully address the question of such a 
cancellation. 

Of course, rapid quenches also displayed interesting (potentially divergent) scal- 
ings in the physical properties of the final state, e.g., the energy density as shown in 
eqs. (5.9) and(5.31). 18 These scalings do not depend on the derivatives of the couplings 
as we set <9(Aa(£ — > ±00) — > to ensure that the system begins at equilibrium and 
eventually settles into a new equilibrium state. Rather these scalings can be associated 
with the behaviour of the expectation value of Oa itself, as follows: Recall the basic 
observation that when the boundary theory has a time- varying coupling Aa(^), the 
stress-energy tensor is not longer conserved, i.e., see eqs. (2.38) and (2.48). Instead, 
the diffeomorphism Ward identity states that for our quenches: 



That is, the time-variation of the coupling is performing work on a system. Of course, 
eq. (6.5) can be written in an integral form as 



Since we only considered quenches where the change in the coupling was finite, the 
(divergent) scaling in this expression for fast quenches must come from {Oa)- 

From the above perspective, the integral appearing, e.g., eq. (3.21) defining a£° 4 
can be seen to describe the total work done on the system. In fact, in the subsequent 
expression for appearing in eq. (3.22), the extra constant terms cancel and so 

we recover precisely eq. (6.6), i.e., an integral form of eq. (6.5). With respect to the 
comments above, recall that the leading scaling in the response was slightly stronger 
for fast quenches than that found for a^° 4 . That is, the integral in eq. (6.6) tempers the 
scaling of the response to fast quenches, as was discussed in detail in section 5. We have 
made some preliminary steps towards extending our present work to quenches with a 
generic relevant (or marginal) operator with dimension A [29]. This work suggests that 
our previous results generalize in a simple fashion with 



18 As noted in section 5, the same scalings appear for the temperature, pressure and entropy density 
of the final state. 



d t £= (O a ) d t \ A . 



(6.5) 




(6.6) 




2 < A < 4, 
A = 2, 



(6.7) 
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for fast quenches with a C 1. Here we might note that the scaling above with A = 4 
matches that found in [13]. 

It is well-known that the renormalization procedure in QFT is scheme-dependent. 
In the context of gauge/gravity duality, such scheme-dependence manifests itself with 
the appearance of finite counterterms in the holographic renormalization. In this pa- 
per, we carefully enumerated all of the finite counterterms for the situation where the 
couplings for O2 and O3 vary in time (or space) for the mass-deformed version of M = 4 
SYM theory. It is worth noting that the scheme-dependence in holographic renormal- 
ization does not introduce any ambiguity in the diffeomorphism Ward identity (6.5) 
- see also eqs. (2.38) and (2.48). That is, in general, both the stress tensor and the 
expectation value (C?a) are scheme-dependent but these ambiguities cancel in eq. (6.5) 
and this Ward identity is still true as written in any scheme. This might be contrasted 
with the conformal Ward identity, where the scheme- dependent terms can make an 
explicit appearance, as shown in eqs. (2.40) and (2.50). 

As our analysis showed, the number of scheme- dependent parameters introduced in 
renormalization of a general deformation (1.2) depends on the choice of A, the dimen- 
sion of the operator. Specifically, we found three ambiguity coefficients in eq. (2.37) for 
A = 3 and two such coefficients in eq. (2.47) for A = 2. 19 We would like to emphasize 
that these ambiguity coefficients are actually necessary to properly interpret the bulk 
gravitational data in the language of the boundary gauge theory. As a simple example, 
consider the case A = 2 with a constant coupling. As illustrated in eq. (2.72), already 
to leading order in A2 in our high temperature expansion, the energy density and pres- 
sure of M = 4 SYM plasma are modified by a term of the form m\ lnT — see [25] 
for more discussion. Further, a similar term appears in the expectation value (0 2 ) as 
shown in eq. (2.75). To make sense of the lnT, one has to introduce an arbitrary scale 
in the theory, as shown in eq. (2.72). The arbitrariness associated with the scale A is 
encoded in the arbitrariness of the scheme-dependent coefficient S\ in eq. (2.47). We 
note that the analogous Si term in eq. (2.37) would produce a similar effect for A = 3, 
e.g., the energy density would be modified by a term proportional to mj lnT, but 
clearly such contributions would only occur at quartic order in the high temperature 
expansion and so were not studied here. 

Of course, the above example illustrates that the time variation of the couplings is 

19 Note that for the case A = 3, this means there are two extra coefficients compared to renormal- 
ization of the theory with static couplings, while one extra coefficient appears for A = 2. 
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not intrinsic to scheme-dependence appearing various natural 'observables'. However, 
once we consider time-dependent couplings, one can expect new contributions to the en- 
ergy density or pressure proportional to In T with coefficients involving time-derivatives 
of the coupling Aa- In particular, the latter must form dimension four operators and 
so can be read off from the first line of eq. (6.4). Hence there are precisely two such 
contributions for A = 3: 

(d t \ 3 ) 2 \nT, A 3 ^A 3 lnT. (6.8) 

As in static case, the proper interpretation of these In T terms requires introduction of 
two new independent scales, A 2 and A 3 , as appear in eqs. (3.20) and (3.31). As shown 
in eq. (3.19), the ambiguity in the choice of these scales is then precisely accounted 
for with the two new scheme-dependent coefficients, 5 2 and <5 3 , appearing in the finite 
counterterms in eq. (2.37). 

Of course, our description of quenches with C 2 contains the same In ^ ambiguity 
which appeared in the equilibrium thermodynamics. However, these quenches also 
contain a novel power law ambiguity arising from the finite counterterm proportional 
to W(p in eq. (2.45). This counterterm contributes a term proportional to the dimen- 
sion four operator (p[ )" ~ cfA 2 in the pressure, as shown in eqs. (3.45) and (3.56). 
However, this term carries a pre-factor Ag (m°) 2 /Tj 4 here. This makes this contribution 
novel in two respects: First of all, the renormalization scale A 2 appears with a power, 
rather than in a logarithm. Note that we are therefore free to set this scale to zero, 
which may be the most natural choice in this case. Second, this term carries a factor of 
(mjj) 2 /Tf and so appears at first order in the high temperature expansion, in contrast 
to all of the other perturbative corrections, e.g., in eq. (3.45) which are second order in 
this expansion. Hence (if A 2 ^ 0), this term formally dominates the behaviour of the 
pressure while the mass coupling A 2 is varying. Here we might note that all of the new 
ambiguities discussed here only play a role while the couplings are varying. Hence, for 
example, they are not important after the coupling has achieved its final value but the 
system is still relaxing to its final equilibrium state. Further, they do not appear in 
expressions relating the properties of the final and initial equilibrium states, e.g., in 
eq. (3.22). 

In standard QFT, the scheme-dependence associated with choosing renormalization 
scales can typically be eliminated with a judicious choice of observables or by choosing 
a physical reference point. For example, calculating the overall vacuum energy density 



61 



in a QFT typically yields a result depending on some cut-off scale and so which is 
scheme- dependent. However, any excited stationary state has a definite energy as 
measured with respect to this vacuum energy. A similar strategy might be used in 
considering the thermal ensemble of the mass-deformed gauge theory with Oi- That 
is, we could choose a particular reference temperature T (7^ 0) and compare, e.g., the 
energy density of a given thermal ensemble to that of the reference ensemble at T . 
This prescription would then yield a result for the equilibrium energy density which is 
independent of the renormalization scale A, 

£(T, T ) £(T) - S(T ) = (V - T 4 + ^ In ^ . (6.9) 

Unfortunately, this approach does not extend naturally to quenches where the mass 
coupling varies in time. However, another strategy that is available, at least in this case 
of the dimension two operator, is to formulate a new scheme-independent observable 
by combining £ and (02). That is, although both of these quantities are individually 
scheme-dependent, as shown in eq. (3.45), the following combination is independent of 
both A 1 and A 2 : 

£{t) ee S{t) + M£ pl ;0 (r) (0 2 (t)>. (6.10) 

In particular, this scheme-independence extends to the situation where the mass cou- 
pling varies in time. However, it is not clear if this definition applies or can be extended 
to higher orders in the high temperature expansion. Further, it is not possible to ex- 
tend this strategy to define scheme-independent pressure. Similarly, no such definitions 
of scheme-independent observables seem possible in the case of quenches with O3. 

Hence it seems like that situation here for the holographic QFT's with time depen- 
dent couplings is similar in many respects to QFT in curved spacetime [27]. In a generic 
background without any special symmetries, one simply finds that (T^) is ambiguous. 
In this context, an 'axiomatic' approach has been developed in which general physically 
reasonable properties are imposed on the renormalized stress-energy tensor, without 
reference to any particular renormalization techniques [27,41,42]. Three of these condi- 
tions would naturally carry over to a general discussion of quenches, or more generally 
time-dependent couplings, in QFT. Namely, the conditions related to causality, finite 
matrix elements for orthogonal states and recovering standard results in Minkowski 
space or here in the limit of time-independent couplings [41]. The first condition, 
however, was covariant conservation of the stress tensor and in the present context 
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Figure 16: (Colour online) The evolution of the a-rescaled normalizable component, 
a 2 px t 2, as a function of r/ot during quenches with the profile given in eq. (6.11) with 
different values of a. 

this would have to be changed to incorporate work terms induced by the time-varying 
couplings, as found in eqs. (2.38) and (2.48). We re-iterate that the latter were dif- 
feomorphism Ward identities, which naturally appear without any scheme-dependence 
in the present holographic renormalization. We might add that the appearance of 
scheme- dependent contributions in the trace anomaly is not resolved by the axiomatic 
approach discussed above [42] and this feature is also evident in the present holographic 
analysis with eqs. (2.40) and (2.50). 

Given the ambiguities in the observables, E , V and (Oa), we demonstrate that some 
care must be taken in defining the excitation and equilibration times. As an example, 
we focus on the equilibration time for quenches with (9 3 and consider a profile of the 
mass coupling which given by 

Pl , 0= Q + itanhl) (l + j^). (8.U) 

With this profile, there is still an exponential turn on of the coupling are early times 
but the approach to the final equilibrium value is now controlled by a power law, i.e., 

2 

Pi q - 1 + + • ■ • as r - > +oo. Fig. 16 shows the response of the normalizable 
component p 12 in these quenches for select values ofa = {l,~,~,|,i}. If we use 
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the same definitions as before, i.e., eqs. (5.13), (5.14) and (5.16), the results for the 
equilibration time scale are essentially unchanged. That is, r eq ~ 1 or t eq ~ l/Tj and 
for r > T" e(J , pi ; 2 follows the adiabatic response (4. 13). 20 Now we may also consider 



defining the equilibration time scale rf directly in terms of a more physical quantity 



■eq 



such as the expectation value of 3 — similar remarks would apply using SorV below. 
That is, we would replace our previous definitions with 

3 (t) - Of iabatic (T) 



^equilibrium 



(6.12) 



where 

equilibrium^ i im a 3 (r) = N 2 T 2 m° f . (6.13) 

r^+oo 2V27T 2 

Then the equilibration time r® q would be defined as the latest time at which 5^ exceeds 
some threshold, i.e., 5®( T > r e ^) < e - Unfortunately, this approach would lead to 
ambiguous results because of the scheme- dependent contribution in the expression for 
3 given in eq. (3.20). In particular, we find 

r0 , s „ . nTi a 2 , TiTi / T Pn \ 4 , 7rT,; 

°ne q ( T ) 



, cx p% In ^ ~ ^ In ^ ~ f ^ j a 2 In , (6.14) 

ambiguity A2 T A2 \ T / A2 



where r eg is the previously defined equilibration time. This ambiguity could substan- 



tially affect the computation of the relaxation time (i.e., making it longer than the 
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thermal time scale) if the renormalization scheme is such that 

a 2 ln^p>l. (6.15) 
A2 

Of course, one may regard the real source of the issue here to be the power law behaviour 
of the profile (6.11) at late times. The latter was introduced precisely to enhance the 
importance of the scheme-dependent terms in eq. (3.20), as these issues would not have 
occurred with the exponential fall-off of the original profile (4.1). Still this example 
stands as a cautionary note in studying the excitation and equilibration times in generic 
time-dependent settings. 

Further it is interesting that, in some cases, the renormalization scheme ambiguities 
can be used to remove the divergent scaling behaviour appears in the one-point corre- 
lators because of the rapid variation of the couplings in these observables. 21 Consider 



20 Wc also explicitly checked that the essentially same QNM behaviour occurs for these quenches 
(6.11) as shown in fig. 8. 

21 These divergences should be distinguished from that in the expectation value of the corresponding 
operator Oa, as discussed above. 
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a quench of (9 3 where the coupling A 3 exhibits a power law growth in the vicinity of 
t > with 

\ 3 (t)=t?, < /3 < 1 . (6.16) 

From eq. (2.36), we can see that for generic values of the scheme-dependent coefficients 
S 2 and # 3 , the observables of the theory would diverge as t — > 0. The divergence 
introduced by the time variation of the coupling p oc A3 is given by 

1 11 

87rG 5 £\ coupling - divergence = - j^Cpo) 2 + 3P0P0 + ^ 2 (p'o) 2 

oct 2 ^- 1 )(3/3-4 + 6/35 2 ), 

oc t 2 ^ _1 )(— 1445 3/ 5 + 725 3 - 3/3 + 2 + 185 2 /3) , 



16 - G 5 ( 3)| coupHn ,_ droer , ence + 2« oc t^(l + 45 2 



Here again, without loss of generality, we have set ax = 0. Hence by choosing 



(6.17) 



12 7-6/3 
52 = -2 + 3^' 53 = 36(2/3 - 1) ' (6 - 18) 

we would remove the divergences in the energy density and the pressure. However, the 
expectation value of the operator O3 would still be divergent. Alternatively, one can 
remove these coupling-induced divergences in the pressure and (O3), with the energy 
density remaining divergent as t — > 0. Presumably, this example illustrates that it is 
more useful to consider the scaling behaviour of observables characterizing the final 
equilibrium state. 

Our analysis revealed that for rapid quenches with TTj 1, the system relaxes, 
almost universally, on thermal scales ~ 1/Tj — see eqs. (5.23) and (5.38). Such a 
response might be an artifact of our treatment of the quenches perturbatively in the 
amplitude of the coupling, effectively about its zero value. It would be interesting to 
study the full nonlinear problem where not only the rate change of the coupling but 
also the amplitude of the coupling (and its variation over the history of a quench) 
can be arbitrary. From figs. 2 and 9, it is clear that even though the relaxation time 
scale for fast quenches of A3 and X 2 couplings is roughly the same, the amplitude of 
the response due to quenches of A3 is quadratically enhanced (in the rate of change) 
compared to that for quenches of X 2 - In this sense, the quenches due to O3 operator 
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more closely follow the intrinsic time-scale of its coupling. Given our perturbative 
analysis of the quenches here, we could treat quenches due to A 2 and A 3 separately. It 
would be interesting to study "supersymmetric" quenches with 

X 2 (t) = X 3 (t) 2 , (6.19) 

which would necessarily require a full nonlinear analysis. 

Another regime of the quenches which exhibited interesting universal behaviour was 
a >> 1, corresponding to slow nearly adiabatic quenches. Recall that for fast quenches 
(a< 1), the normalizable mode of the bulk scalar relaxes to its final equilibrium with 
a quasinormal mode behaviour, e.g., see figs. 8 and 15. In contrast, the relaxation 
in the nearly adiabatic regime occurs without excitation of any quasinormal modes. 
Rather for a > 1, the deviation from the purely adiabatic evolution in eq. (4.13) 
is governed entirely by derivatives of the source profile, e.g., a d T pifi(r/a) for O s 
quenches or by a d T p l 10 (r/a) for 2 quenches. Indeed, in either case with an arbitrary 
source p = p (r/a) and a >> 1, we can solve the bulk scalar wave equation (3.3) 
perturbatively in 1/a: 

oo 

0! =po(r/a) 4>\(p) + a ~ n tfV/a, p) , (6-20) 

n=l 

where 4>\{p) is the equilibrium solution normalized such that the leading coefficient in 
the asymptotic expansion is one. For example, for the m 2 = —3 scalar, the latter would 
be the solution given in eq. (2.55) divided by a factor of A/i. Now for n — 1, we find 

<pf\r/a,p) = ad T p (r/a) g^(p), (6.21) 
where g^ satisfies the inhomogeneous equation 

4(1 " P 4 ) *„9 W + "-if V 1} + $ 9 (1) = ^l~ Ml ■ (6-22) 

Similarly, successively higher derivatives of p appear at higher orders in this expansion. 
Hence the entire dynamics of the linearized solution is simply governed by the nearly 
adiabatic driving source. 

Quantum quenches from a thermal initial state were discussed previously in [5]. 
The major difference between this seminal study and our present analysis is that we 
are considering a strongly coupled field theory while the authors of [5] investigated 
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quantum quenches of a free scalar field. Effectively, the momentum modes in the latter 
theory reduce to a collection of decoupled simple harmonic oscillators (SHO's). In [5], 
it was shown that a(n instantaneous) quench of an individual SHO from an initial 
frequency Wj and temperature Tj to a final frequency Uf results in an effective thermal 
final state with temperature Tf 

T f = \ [\ + ^ Ti , for Ti > . (6.23) 

A striking result here is that if we also choose Uf <C Wj, then the final temperature 
is one-half the initial temperature. Quite generally ujf < Wj results in Tf < Tj, i.e., 
the system cools! A similar cooling effect was found for quenches of the mass in the 
field theory where Tj ^> m 8 ^> mf in , using an effective temperature determined by 
averaging over all of the momentum modes. In fact, Tf = T/2 was found for such 
quenches in d = 2 and 3. However, in d — 4, the effective temperature was dominated 
by high momentum modes and they found Tf = Tj. 

Recall that our analysis also applies to the same high temperature regime, e.g., 
T ^> m°f b . Hence it is interesting to compare the above results to those for our quenches 
from the massive gauge theory to the CFT, in which case we also satisfy raj ^> mf in . 
The expressions for Tj/T are given in eqs. (3.33) and (3.58) for these quenches with 
3 and 2 , respectively. Recall that in section 3, we argued that a^ 4 ~ for a ^> 1 
- see the discussion around eqs. (3.26) and (3.51) — and hence it is clear that these 
quenches are indeed cooling the system provided the transition was sufficiently slow. 
However, as the characteristic time scale of a quench becomes shorter than a typical 
thermal scale 1/T, we also saw that af A grows large and so the corresponding quenches 
always result in an increase of the temperature. The precise thresholds between the 
cooling and heating regimes are given in eqs. (5.6) and (5.29). Of course, since our 
analysis applies for the high temperature regime, any changes in the temperature are 
perturbatively small in m°j/Tj and so much like the result of [5] for d = 4, we have 
Tf ~ Tj to leading order. Hence we may conclude that the thermal quenches studied 
in [5] for a free field and those studied here for a strongly interacting QFT share certain 
similar features. Of course, we pointed out above that for fast quenches, the behaviour 
of these two systems is qualitatively dissimilar. 

We would like to contrast our framework with previous analytic investigations using 
the AdS-Vaidya solution [10]. These studies investigated issues related to thermaliza- 
tion of a holographic plasma by considering the gravitational collapse of a thin shell 
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of 'null dust'. Of course, the latter scenarios might also be regarded as a holographic 
quantum quenches. However, we would note that by design, the null dust does not 
excite any local gravitational degrees of freedom, in particular, on the horizon. For 
example, the spacetime becomes that of a static black hole immediately beyond the 
boundary where the bulk stress tensor goes to zero. In the dual description then, this 
scenario injects energy into the gauge theory with an 'exotic' probe which does not ex- 
cite any of the local degrees of freedom in the holographic plasma. Hence this approach 
may be useful to study questions related to the formation and growth of the horizon, 
however, it can not address many aspects of the holographic quench. For example, it 
does not reveal the scaling behaviour of the response or the quasinormal response of 
the horizon. Our analysis indicates that both of these features play an important role 
in the equilibration of the system. It is interesting that the Vaidya solutions seem to 
be a good approximation to a quench induced with a massless bulk scalar in a small 
amplitude expansion [9], similar to that used in our analysis. The analytic approach 
of [9] also allows one to systematically improve the gravitational description within this 
expansion. 

There are numerous aspects of the holographic quenches described here that are 
left for future analysis. First and foremost, we would like to extend our perturbative 
treatment to a full nonlinear analysis of the gravitational equations of motion. How- 
ever, let us outline a few other potentially interesting extensions: 

■ In this paper, we only considered the response as measured by a few simple one-point 
correlators. These are representatives of local physical observables in the holographic 
theory. An open question remains to analyze the behaviour, in particular the re- 
laxation, of non-local physical observables, such as higher-point correlation functions, 
Wilson lines or entanglement entropy. Previous studies [10] indicate that these non- 
local observables equilibriate more slowly than the local observables and so it will be 
interesting to see if this effect persists in the present framework which accounts the full 
quasinormal response of the horizon. 

■ Our present analysis was restricted to quenches which are spatially isotropic and 
homogeneous. As a result, we eliminated 22 some of the potentially important relax- 
ation channels of the system, namely, the sound and the shear modes in off-equilibrium 
gauge theory plasma. In a more general anisotropic or inhomogeneous quench where 
such modes are excited, one may expect that the hydrodynamic modes will have much 

22 We would like to thank Rob Leigh for raising this issue. 
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slower decay rates than the quasinormal modes of the bulk scalar found in the present 
work. Inevitably then, these new modes would dominate the late-time behaviour of 
the system, however, their influence on the equilibration time would still depend on 
how efficiently they are produced during the quench. We defer the consideration such 
inhomogeneous quenches to future work. 

■ We focused on the properties of the energy density and the pressure of gauge theory 
plasma undergoing a quench. We pointed out that both of these observables are renor- 
malization scheme dependent. It is interesting to address the question of ambiguities 
in the entropy density away from equilibrium [43]. In particular then, it would be 
interesting to explore how a choice of a renormalization scheme affects the definition 
of the entropy density. 
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A Discretization of the evolution equation 

We describe here the basic strategy employed to discretize the evolution equation. The 
schematic structure of such equation, for a field /(r, p), takes the form: 



for suitably defined multiplicative factors A..., and a source term Jo, all explicit func- 
tions of (r, p). Notice the principal part of the system, i.e., 



is a second order partial differential equation admitting two propagating modes. One 
of them has a characteristic speed given by A pp while the other propagates instantly 



= d Tp f - A pp d 2 pp f - A T d T f + A p d p f + A f + J , 



(A.I) 



= (d Tp - A pp d 2 pp ) f , 



(A.2) 
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inwards from the boundary (a consequence of having chosen incoming characteristics 
to foliate the spacetime). Further, notice that for a scenario with a black hole present 
from the start, as is the case for our studies here, A pp > (< 0) outside (inside) the 
horizon. Thus only one boundary condition can be prescribed at the AdS boundary. 

A natural strategy to discretize the system, taking advantage of the characteristic 
structure explicitly displayed by equation (A.l), has been presented in [44] (and further 
extended in [45] to include dissipation and in [46] to allow for adaptive meshes). In our 
implementation, we follow this approach to obtain a second order accurate implemen- 
tation in the following way. We introduce a uniform discrete spatial grid covering the 
domain p G [0, L p ] with N p + 1 points located at pi = (i — l)dp (dp = L p /(N p ). Dis- 
crete, evenly spaced, time levels are introduced as r n = r Q + (n— l)dr (with dr = adp). 
The continuous function f(r,p) is then represented by values on this grid structure as 
fij ni p,i) = Only two time levels are required for discretizing our target equation 
as follows. Consider / is known at level n, i.e., all values /" are known, as well as 
those at level n + 1 for p < pi, to obtain the value of a centered discretization at 
the virtual point (n + 1/2, i — 1/2) is achieved by discretizing the terms involved as, 

fn+l fn+l fn _i_ fn 

a f.in+1/2 Ji Ji—1 Ji ' Ji—1 ;i n\ 

drpfli ~ dpTr ' {A - 6) 



fn+l £■ n+1 , rn+1 \ fn o f n _i_ f n 

° p j\i - 2dp2 , {*-V 



rn+l _|_ £n+l en £n 

p, j.in+1/2 Ji 'Ji—1 Ji Ji—1 /\ r\ 

d T }\i ~ ^ , (A.5) 

fn+l fn+l I fn fn 

Pl r\n+l/2 _ Ji Ji-1 ~t~ Jj Jj-i ( \ a\ 

Cp/li " ^ ' { } 

fn+l _|_ fn 

fl n+l/2 = h-1 + h (A?) 

The coefficients A... and the source J are evaluated analytically at the virtual point. 
Plugging these expressions into the equation and solving for defines the scheme 
employed. Notice the resulting scheme is thus an inwards "marching" algorithm that 
provides the value of Further the "slanted" nature of this scheme (i.e., the 

forward offset of the second derivative at level n) is simply a reflection of the charac- 
teristic structure of the system. Since initial data is given at an initial hypersurface, 
and boundary conditions at p\ — 0, this scheme suffices to integrate all points with the 
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exception of p 2 , Pn„+i where the simple change, 



fn o fn I f n 

o2 f in+l/2 _ J3 Z J2 -r Jl i a , 

C,/| 2 - ^ , (A.i 



f n+l f n+1 _|_ ^n+1 \ fn _ o f n _i_ f n 



suffices to provide a first accurate order approximation which does not affect the overall 
second order convergence of the solution. 

We specify the characteristic initial data f mitml (p) a t initial time r = t\. 

f(r = r 1 ,p) = r Mal (p), (A.fO) 

and the boundary data f boundar v a t p\ =0, 

f(r,p = p l ) = f boundar y(r). (A.ff) 

As a representative test of the code's convergence behaviour, we consider the bosonic 
case and obtain solutions for different grid sizes N p = 2 P x 100 (p = 0,...,4) and 
simulate the system until equilibration. We monitor the (L 2 norm of the) difference 
of solutions with successive values of p (defining e p = ||0 p+ i — (f> p \\2) as wen as the 
convergence rate q calculated as, 

2q = ||0 P +i -0 P ||2 = _^v_ ( A 12 ) 

||0p+2 - 0p+l||2 e p+ i 

The results are shown in figure 17, the left panel illustrates how e p decreases as the 
resolution is improved while the right one the convergence rate q illustrating second 
order convergence behaviour. 

As a last remark, we stress that in the general case metric variables governed by 
an evolution equation (e.g., £ determined by eq. (2.14)) can be integrated using the 
same strategy described above. Thus, this algorithm together with straightforward 
integration along the radial direction (e.g., for A determined by eq. (2.15)) suffice for 
implementing the complete system provided consistent boundary conditions are pro- 
vided at the boundary. For other strategies see [13] (employing also a characteristic 
formulation of the equations) and [16,18]) (for works implementing a Cauchy formu- 
lation). 
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Figure 17: (Colour online) Convergence tests. (Left panel) L2 norm of the difference 
between solutions obtained halving the discretization length. (Right panel) Estimated 
convergence rate from the numerical solution, indicating the numerical solution indeed 
converges to second order as expected. 
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